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Chapter 1 



Overview 

1.1 Conventions 

First, a basic note on convention. For what follows, we will be taking the signa- 
ture (— , +, +, +) for the metric tensor and the convention (V a Vf> — V&V a )f c = 
R a b c d Vd for the Riemann tensor. We will also take e i23 = +1- If a mathemat- 
ical quantity (i.e., R ab ,V a , etc.) is written unmodified, it will be considered to 
be a 4-dimensional object. If the same object has a tilde placed over it (i.e., 
Vj), then that object will be considered to be an object intrinsic to a null 
hypersurface embedded in a Minkowskian 4-space. If that same object has an 
overbar (Rij), then that object will be considered to be an object intrinsic to 
a 3-dimensional spacelike submanifold. If the object in question is denoted 
with an overhat (Rab), then that object will be considered to be intrinsic to 
a two-dimensional submanifold to the 4-dimensional space. We will also take 
the convention that indices that run over 4-dimensional labels will be labeled 
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with lower case Latin letters from the beginning of the alphabet (a,b,c,...). 
3-dimensional indices labeled in the intrinsic coordinates of the 3-dimensional 
space will be labeled with lower case Latin letters running from the middle of 
the alphabet (i,j,k,...). Finally, 2-dimensional indices running over the coordi- 
nates intrinsic to a 2-dimensional space will be denoted by capital Latin letters 
running from the beginning of the alphabet (A,B,C,...). Note that, with the 
use of the inclusion operators to be defined below, it is possible to have (for 
instance) a 2-dimensional quantity whose indices are 4-dimensional. 



1.2 Scope of Dissertation 

Over the past three and a half decades, there has been a broad array of work on 
formalizing and redefining the notion of the black hole horizons, starting with 
the work of Regge and Tietelboim [ ]. The primary thrust of this work has 
been to take the heavy mathematical machinery of immersions and manifolds 
and to use them to obtain a broad swath of powerful results in Relativity. 
However, as this has progressed, there has only recently been a focus upon 
how to compute and analyze these quantities. And while there are numerous 
attempts to use these quantities to set boundary conditions for numerical 
simulations [ | j | [4] [5] [6] and for quantum gravitational computations [7], 
there has not been a clear method to use 3+1 conditions in order to track 
the motion of a dynamical horizon across a 3+1 slice. In section 8.1.1, this is 
precisely what is done. It is also explicitly shown, using the Vaiyda solution, 
how this term can be used to track the area growth of a black hole apparent 
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horizon without assuming any prior knowledge about the dynamics of the black 
hole. 

This work, therefore, will have two goals. The first one will be to re- 
derive existing results in a way that makes it transparent how one can do 
concrete calculations using the listed techniques. Theorems from differential 
topology and geometry will be used with a focus towards indicating how one 
can choose a particular coordinate system and calculate quantities. In par- 
ticular, with null geometry, this question is treated with great care, as a null 
surface will have a tangent space and a cotangent space that are not metrically 
related to each other. Thus, this work has many explicit examples and con- 
crete calculations, both in the body of the paper and in the appendicesf. The 
explicit null decomposition of Minkowski spacetime and of the Kerr spacetime 
done in Appendix A is, to the author's knowledge, not currently published 
in the literature. The formalism used in doing this decomposition actually 
improves existing literature, as it resolves an ambiguity in the definition of the 
degenerate metric tensor and the metric compatible null connection noted by 
Ashtekar et al. [8] 

Additional original work provided here includes the analysis of the 
McVittie spacetime 1 . Discussion of the effects of the expansion of the uni- 
verse on geodesies is described, as well as actual numerical computation of the 

geodesies. There is a discussion regarding the singularity of this spacetime 

1 The McVittie spacetime fuses properties of the Schwarzschild and Robertson- Walker 
spacetimes, and thus can be considered to describe a black hole sitting in a Robertson- 
Walker background. 
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at the surface r = 2 M, where curvature invariants C a b c d,C , R a bR and R 
are computed in order to establish the still-disputed claim that the McVittie 
"horizon" represents a spacetime singularity, and not a coordinate singularity a 
la Schwarzschild. Then the spacetime is used to generate results related to the 
Hamiltonian dynamics of Relativity, in particular, as an example of how one 
may make sense of an ADM mass even in a non-asymptotically flat spacetime. 

This work also treats the elegant classical counterterm treatment of 
the ADM quantities, as originally computed by Hawking et al. rather than 
the more standard treatment which rejects boundary terms in a variational 
principle, and then inserts them at the last step in order to ensure consistency 
of the Hamiltonian equations of motion with the Hamiltonian [ ]. While the 
treatment of Hamiltonian Relativity using counterterms is not original, the 
discussion to follow does provide a new spin upon the discussion, focused 
on the interaction between boundary degrees of freedom and bulk degrees of 
freedom. 

Finally, this work derives boundary conditions upon the lapse function 
and a component of the shift vector, and does so both in the intrinsic language 
of the Isolated and Dynamical horizons of Ashtekar and in the 3+1 Hamilto- 
nian language that is best suited for numerical simulations of Einstein's equa- 
tion. And this is all done with an aim of developing a clear technique to show 
how the relevant objects can actually be evaluated and computed in a 3+1 
context. 
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Chapter 2 

Basic Definitions and 
Decomposition Process 

2.1 Choice of a 3+1 slice and normal vector 

In the literature, there has been a dearth of clarity on the topic of projecting 
tensors from higher dimensional spaces onto lower dimensional spaces. In this 
section, a consistent procedure for doing this will be described, and then a 
direct application to the derivation of the (3+1) ADM split will be given. The 
key point to remember is that the projection operation has different behavior 
depending on whether one is using the projection map or the inclusion map 
and whether or not these maps are acting on vectors or one-forms. 

Now, consider a decomposition of a general manifold M with metric 
tensor g a ^. Let M have the local topology Ixm, where m is a submanifold to 
M. Now, let r be some function defined on M which takes a constant value on 
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each fibre m of R. Then, V a r will be a normal vector to the fiber. For now, 
assume that V a TV a r 7^ 0, (a case that will be dealt with on page 72), and 
define £ as the sign of V a rV a r. After making the definition a = ^===^, it 
is simple enough to define the unit normal to m by the formula 

n a = aV a r. (2.1.1) 

Typically, the next step would be to go and define the normal projection 
operator 7^ = g a ^ — £,n a nb, which we will soon do, but some care should be 
taken with this definition, as there are two different senses in which we can take 
the phrase "3-metric". First, we wish to choose a special coordinate system for 
our manifold. In particular, we wish to choose r as one of the coordinates in 
the neighborhood of one of the fibers. Clearly, this should be a valid coordinate 
choice, as r is constant on each fiber, so advancing in r merely advances from 
one fibre to the next, just as stepping forward in t in standard Minkowski 
spacetime merely advances one from one spacelike copy of M. 3 to the next. 
So, after replacing one of the original manifold coordinates with r we have 
a coordinate system of the form (r, x l ), where the index i ranges over the 
remaining coordinates in the spacetime. 1 Therefore, it is now clear that, by 
assumption, we have V a r = (£, 0, 0, 0) 2 , and a = ^^w- Therefore, this 
gives a relatively simple coordinate expression for the raised form of n a : 

1 While the language of this example is geared for the most common case, a timelike V a r 
in a 4-dimcnsional Minkowskian spacetime, and the argument will advance as such, note 
that no such assumption is being made here. The above methodology will work for any 
non-null V a T. 

2 The factor of £ is inserted so as to make g ab W a T future pointing 
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We therefore have an expression for the raised operator 7 in terms of 
the original metric components: 



Where we use the fact that £ = ±1 to set £ 2 = 1. It can now clearly be 
seen that all of the components of / ~f ab where either index takes on the value r 
vanishes. It is therefore natural to consider ^y ab to be an operator that takes 
four-component one-forms on M and projects them onto three-component vec- 
tors on m. Therefore, we can consider the operator r y ab in two senses: First, we 
can think of it as a projection operator living in a 4-dimensional space whose 
output is vectors with three nonzero entries. Secondly, we can think of it as 
an intrinsic 3-dimensional object that serves to define the inner product over 
the covector space to m. Now, we will show that the natural inverse of the 
three by three matrix 7^ is in fact obtained simply by restricting the 4-metric 
to the appropriate 3-dimensional indices, g^: 




(2.1.3) 



f 3 9jk = (9 



1 



)9jk , 



(2.1.4) 
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• TT 



1 



9 Ti (g Ta 9ak - g TT g T k) , 



(2.1.5) 



= 9 ' 9ak~ 9 ' 9rk~ 



• TT 



9 



S\ - 9 %T 9rk - — 5 T k + g Tl g Tk 



9' 



5\.. 



(2.1.6) 
(2.1.7) 



Here, we change the summation over 3-dimensional indices to a summa- 
tion over 4-dimensional indices by the identity v l V{ = v a v a — v T v T , and, in the 
transition to the last line, we realize that k ranges over only spatial indices, so 
therefore 5 T k is identically zero, and we use the symmetry of the metric tensor 
to cancel the first and third terms. Therefore, we can see that the projection 
operator on covariant indices merely requires that we drop the r component of 
the one-form, and this is where care should be taken, as if we use g ab to raise 
and lower the indices of / y ab , we could na'ively believe that the proper three 
projection of a one-form uj a might be 7 a (,c<; a 3 . A simple computation, however, 
will show that this form does, in fact, have a r component, and therefore, is 
not properly a 3-dimensional one-form, in the sense that it is three quantities 
along the three coordinate basis one-forms on the cotangent space to m. In- 
stead, as the above example involving the metric indicates, the proper thing 
to do is, instead, simply drop the dr component of the covector u a . 

Finally, we might want to project a 4-dimensional vector onto our 3- 
dimensional space. This action, however, has a more complicated form than 
simply dropping the r component. In order to accomplish this, we first lower 
the index of the vector, follow the above projection procedure, and then raise 

the vector with 7 y : 

3 And, as we will see on page 12, in a certain sense it is 
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v a ^g ia v\ (2.1.8) 

g ia v a =g iT v T + gijV j , (2.1.9) 

v k =-y ki girV T + T ki g ij v j , (2.1.10) 

=(^ + nV) ffitU t + 7j -V, (2.1.11) 

= (5 T * + n fc n> T + 7j V , (2.1.12) 

=7aV. (2.1.13) 



So, while the operator 7a 6 does not take four component one forms on 
M and spit out three component one forms on m, it does do this to vectors. 

On the other hand, it should be clear that we might have a problem 
going in the other direction, from 3-dimensional space to 4-dimensional space. 
In particular, if we were to map from the 3-space back into the 4-space, we 
would want our target vector/one- form to be normal to the appropriate version 
of n a , since we want to consider n a to be the normal to the 3-space. But, if we 
were to merely map a general 3-dimensional one form in the trivial manner: 

Si = (x 1 ,x 2 ...,x n ) ->■ s a = (0,x 1 ,x 2 ...,x n ) . (2.1.14) 

We can see right away that we will get a result that is inconsistent with 
our desire to be normal to n a , since n a = ~^— g ta and, in general, the g t% ^ 0. 
So, what are we to do? The answer is that, for one-forms, we use the operator 
7 a fe not as a projection operator from M onto m, but rather as an inclusion 
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operator of m into M. If we use the above mapping, and then follow that with 
having the one-form be acted upon by 7%, we get a vector that is normal to 
n a , since 

n a la b = n a (5 a b + n a n b ) = n b - n b = . (2.1.15) 
Furthermore, we have 

l ab lbc = (g ab + n a n b )( lbc + n b n c ) = 6 a c + n a n c + n a n c -n a n c = 5 a c + n a n c = 7 a c . 

(2.1.16) 

Now, we will show that use of 7 ab enables us to find a vector in M that 
is, in a sense, equivalent to a one form living in the cotangent space of m, 
thereby avoiding this problem. Then, once we know this duality, we can then 
use the appropriate metric tensor to raise and lower indices as we will. First, 
consider a one-form Sj living in m. Define s a according to the trivial inclusion 
map seen in (2.1.14). Since we know that the components of are equivalent 
to the components of 7^, and that the matrix ^ is the inverse of the matrix 
jij, while, at the same time, we know that 7a b satisfies the identity 

l ab l bc = {g ah -in a n b ) {g bc -in b n c ) , 
=5 a c - in a n c - £n a n c + ^n a n c , 
=5 a c -£n a n c , 
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=7V 



(2.1.17) 



Then, we can show that if we define s a = 7° b s b , we can show that the 
norm of S{ under 7^ is equivalent to the norm of s a under g ab : 



s^^SiSj., (2.1.18) 

=l ab s a s b , (2.1.19) 

=lc d l ca l db SaS b i (2.1.20) 

= (^-^c^)7 ca s a 7 d6 s & , (2.1.21) 

=9cdS c s d , (2.1.22) 

=s a s a . (2.1.23) 



So, then, we can lower the left hand index of ^ ab using g ab , and then, we 
can use the result, 7^, as an inclusion map whose target has the same norm in 
M as the source had in m. Also, note that since n a has only a r component, 
we can simply drop the r component of s a to recover Sj. Henceforward, this 
paper will drop the bar notation on one forms, and let the type of index denote 
whether we are discussing the 3-dimensional one-form or its inclusion into the 
four-manifold. 

Meanwhile, using the trivial inclusion map for vectors creates no prob- 
lem, since V a T and n a have no components along any other direction than the 
r coordinate. Therefore, if we take: 
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s i = (x\x 2 ...x n ) ^s a = (0,x\x 2 ...x n ) . (2.1.24) 

We are guaranteed to not run into any trouble, since the right hand of 
the above equation is normal to all of the one-forms with no components that 
project along the space, and it is guaranteed to have the same norm, since 
9ij = lij- So, we have the following rules for the projection and inclusion 
maps: 

1. If we are using the projection map: 

(a) One forms project trivially: just drop the appropriate component, 
and the projection is done 

(b) Vectors, however, do NOT project trivially, and need to be pro- 
jected using the 7% operator, which is equivalent to lowering with 
the full metric, operating with the trivial projection map, and then 
raising with the three-metric 

2. If we are acting with the inclusion map (we are taking 3-component 
objects intrinsic to a 3-surface, and seeing how they look in the full 
spacetime) 

(a) in this case, it is the vectors who have the trivial behaviour. Just 
add a zero in the appropriate component, and then you have a 
vector in the full tangent space 
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(b) One forms, however, need to have the inclusion matrix / y a b oper- 
ate on them. This is ultimately equivalent to raising them using 
the induced metric j %3 , operating with the trivial inclusion map 
for vectors and then lowering with the enveloping four-dimensional 
metric. 

The above procedure is justified by a look at any differential topology 
textbook (see, e.g., [9]), which will define the action of the pullback operator 
(that maps from a manifold with higher dimension to one with lower dimen- 
sion) in terms of the differential of the map from the higher dimensional space 
to the lower dimensional space. It will therefore have a trivial action on forms. 
Meanwhile, the book will define the inclusion operator (which acts in the oppo- 
site direction) in terms of the vectors in the target space that are being acted 
on by the above forms, and therefore, the inclusion operator will act trivially 
on vectors, but not forms. The topology book will typically stop short of 
saying what to do with the non-trivial actions, as a metric tensor is typically 
not assumed. In our case, however, we can use the metric and induced metric 
to move the tangent space into the cotangent space, and then operate with 
the inclusion map and the projection map, and then move back between the 
tangent space and cotangent space to get the correct operator. 

So, in summary, there are two different senses in which we can talk about 
3-dimensional tensor indices. The first is to consider the tensors as arising 
from three component indices intrinsic to the 3-dimensional submanifold. The 
other is to consider the tensors as being the targets of the inclusion map 
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described above. These objects have four indices, and live in the 4-dimensional 
spacetime, but are "3-dimensional" in the sense that they are normal to the 
vector n a and also to the one-form n a . Above, we have shown that the two 
methods are equivalent, and have shown the manner in which we can map from 
one version of a 'three- vector' to another. From here on out, we will work 
primarily with the version of a three vector that lives in the 4-dimensional 
space, and then project down when necessary. 

Finally, one might object to our above procedure, since we took the 
ten metric tensor components, and projected them down to get the lapse a, 
as defined in 2.1 and the three-metric jab as defined in (2.1.3), which has 
six independent components. What of the other three 4-metric components? 
These three components are encoded in the shift vector, defined according to: 



It should be clear that this vector does, in fact, have three components 
that correspond to the metric tensor components that are dropped when the 
projection is done. And with that, we have encoded all of the ten independent 
metric tensor components in terms of equivalent 3-dimensional quantities. At 
times, it is convenient to speak of the time evolution vector, which is the 
vector that generates the evolution in the three-geometry as one advances 
along surfaces with increasing values of r: 



(2.1.25) 



t a = an a + 7 a % . 



(2.1.26) 
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In fact, it is easy to show that, in our coordinate system, this vector 
has only a r component: 



t a =an a + 7 a6 /3 6 , (2.1.27) 
+ 7 a V, (2.1.28) 



1 N ^ „ab. 



fa 7 

J +g ab g b r-Zn a n b g bT , (2.1.29) 



TT 

ra 



- g — + 5\-in a n T , (2.1.30) 

ra £ to. t 

= ^- + 5 a r -^(4^=)^^, (2.1.31) 

--5\ + - y?—, (2.1.32) 

g rr £gTT 

--5\. (2.1.33) 



Now, let us consider the appropriate way to define a 3-dimensional 
derivative operator. We remember from basic differential geometry that a 
derivative operator is unique if it (i) has domain and range within the manifold, 
(ii) maps the metric tensor to zero, and is (iii) torsion free. Therefore, consider 
the following definition, for any vector v a such that v a n a = 0: 

VaV^Yalt^cVd. (2.1.34) 

Condition (iii) is automatically satisfied if the original connection is tor- 
sion free, since you can't introduce antisymmetry in the Christoffel symbols by 
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multiplying them by matrices. Condition (i) is satisfied due to a combination 
of the restriction on v a and equation (2.1.15) above. Appealing to the metric 
compatibility of the original connection, we can explicitly show that condition 
(ii) is satisfied: 

V a 7fe c =i m ai n bl k cVm(gnk + n n H k ) , (2.1.35) 

=0 + l m al n bl k c{n n W m n k + n k W m n n ) , (2.1.36) 
=0. (2.1.37) 

Since 7% annihilates n a . Therefore, V a is the unique derivative operator 
on m. Furthermore, we can note that this is, in fact, just equivalent to the 
standard definition one would make of a covariant derivative, if we first note 
that the operator ^ a b merely serves to project out any derivatives along the 
t direction. We can therefore take r ji a 'jj b d a Vb = d(Vj after having dropped 
the t indices according to the procedure outlined above. Now, we can check 
the equivalence of the definition (2.1.34) with the standard definition of the 
covariant derivative: 

7* VVa^ =7* V(«^ - T ab c v c ) , 

=9 i V i ~ ]pi a lj b Vkl k c9 Cm {9am,b + 9bm,a ~ 9ab,m) , 

=diVj - \li a lj b v k '~i km [d a ('~i bm - n b n m ) + d b (j am - n a n m ) 
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d k (lab -n a n b )\ , 

3 2 [iam,b ~r lbm,a iab,m\ i 



=diVj - Yij k v k 



Where, in the transition from the third to fourth lines, we saw, that 
when the partial derivative operator acted on the n m n n terms and we expanded 
with the product rule, we saw that it was impossible to not contract one of the 
n m on a corresponding 7 m „, yielding zero. Similarly, when going to the last 
line, we realized that, if all of the quantities involved are 3-dimensional, then 
the 7 a 6 operator is merely the identiy operator on the 3-dimensional space, 
and we defined the f ab c operator to be the obvious Christoffel symbol of the 
3-dimensional space. 

Now, beyond just the basic metric and connection of the three manifold, 
we should consider the information about the embedding of m in M. In 
order to do this, we introduce the extrinsic curvature tensor according to the 
definition: 



Where the second factor of 7^ can be dropped in favor of 5f since those 
two tensors differ only by a factor of n&n d , which vanishes when contracted on 
the covariant derivative term since n a has a fixed norm. This concludes our 
basic discussion of how to encode 4-geometry in terms of 3-geometry. In the 
next chapter, we will investigate how to translate the curvature and Einstein's 



K ab = -lilt y c n d . 



(2.1.38) 
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equation into three-dimensional language, using a more explicit (albeit less 
elegant) methodology than is typically used in the literature. 
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Chapter 3 

ADM decomposition of Einstein's 

Equation 

In this section, we will now proceed to take the Einstein Equation, and express 
it in terms of the 3-dimensional quantities described above. This can be (and 
in the literature, usually is) done by writing the Hilbert action f Rd 4 x 
in terms of the 3-dimensional quantities, defining the conjugate momentum 
to '-fab, taking a Legendre transformation of the Lagrangian in order to find 
the Hamiltonian, and then writing down Hamilton's equations of motion. In 
this chapter, we take a different approach, explicitly calculating Lie deriva- 
tives of appropriate quantities using a 3+1 split. While this approach is 
much less elegant than the Hamiltonian technique, it serves a few purposes. 
First, it gives desired results much more directly-there is no need to Legen- 
dre transform, define conjugate momenta, and the like. In the Hamiltonian 
approach, the equations you naturally get involve the canonical momentum 
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U ab = (pi ab K — K ab ) , and the equation for IP 6 must be used to derive 
an equation for K ab , which involves significant amounts of algebra involving 
factors of A f a j ) . In the brute force method, with little embellishment, one can 
directly calculate 7 a b and K^. Second, while both techniques are available 
for splitting a Lorentzian 4-manifold of local topology M x m, only the brute 
force technique used here is generalizable to the splitting a neighborhood of 
the boundary of a Riemannian 3-surface into its 2-dimensional boundary and 
a 1-dimensional normal space. 

For completeness, and for the sake of investigating the boundary terms 
that the Hilbert action necessarily gives rise to, we will compute the ADM 
equations using the Hamiltonian picture in Chapter (4) 

3.1 Constraint Equations 

First, given an arbitrary four vector v a , let us recall the definition of the 
Riemann Curvature tensor: 

R abc d v d = 2V [a V b] v c . (3.1.1) 

Since we defined the 3-dimensional connection above in equation (2.1.34), 
we also have a simple enough expression for the 3-dimensional curvature tensor, 
of an arbitrary 3-dimensional one-form u c (and using the metric compatibility 
of the 4-connection, combined with n a ^ b a = 0): 
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R abc d u d =2V [a V b] u c , (3.1.2) 

=2 1 m [al n b]1 k c \/ m (l r nl S k^rCo s ) , (3.1.3) 
=27 m [a 7 r 6] 7 s c V m V r a; s + 27 m [a 7 r fe] 7 A: c (V m 7A : s )V r w s 

+ 2 1 m [al n h]1 S c{V mln r )V r u s , (3.1.4) 
=la m lb n lc k R m nk d uj d - 2£ 1 m [al r b]l k cn s (y m n k )V r uj s 

- 2^ m [al n b]1 s c n r {V m n n )V r uj s . (3.1.5) 

Before we proceed farther, we should derive some properties of the gra- 
dient of n a . First, as asserted above, we know that n b V a n b = 0, since it is easy 
to show that it is equal to minus itself. Now, let us evaluate n a V a n b . Using 
Equation (2.1.1): 



n a V a n b =n a V a ^aV b r) , (3.1.6) 

=in a V b rV a a + £an a V a V b T , (3.1.7) 

=n a — V a a + £an a V b V a T, (3.1.8) 
a 

ii 

=n a n b VJn(a) + an a V b (—) , (3.1.9) 

a 

=n a n b V a ln(a) - oin a n^-^- + n a V b n a , (3.1.10) 

cr 

=n a n b VJn(a) - £V b ln(a) , (3.1.11) 

= - £V 6 Zn(a) (5 b a - £n a n b ) , (3.1.12) 
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= - il bc l ca V a ln{a) , 



(3.1.13) 
(3.1.14) 



Finally, we can see that, since 7 annihilates all vectors parallel to n a , 
we have 



7 am 7 6n V m n n = 7 am 7 bn V m («V n r) = a 7 am 7 6n V m V n r . (3.1.15) 
By the vanishing of the torsion tensor, K ab is symmetric 1 , and we have: 

V a n fe = £n a V fe In (a) - K ab . (3.1.16) 



Using these results, we can see that the second line of equation (3.1.5) 
vanishes, since it is a 3-dimensional antisymmetrization of (3.1.16). Using 
the fact that n a u a = 0, we use the product rule to replace (3.1.5) with the 
following. We then proceed by factoring out the factor of u s and multiplying 
by isd- 



Rabc d CO d =^7cRmnk d C0 d + ^TlUc^s(V m U k )V r Tl S 

~ ilTl r a l k c uj s (y m n k )V r n s , (3.1.17) 

1 Note that although K ab is symmetric, V a nb clearly is not, as should be explicitly clear 
in equation (3.1.16). It is only the projection of V a n;, onto the 3-surface that is symmetric. 
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(3.1.18) 



Rabcd =la m lb n lc k l l dRmnkl ~ £K ad K bc + ^K ac K bd , 



(3.1.19) 



R = 7 ac 7 bd Ra bcd ~ iK ab K ab + iK 2 , 



(3.1.20) 



R - 2£R ab n a n b - iK ab K ab + iK 2 . 



(3.1.21) 



Having come this far, now we have to assume that we are making a 
timelike slicing of our spacetime. For the rest of this section, let us therefore 
assume that £ = — 1, meaning that n a is a timelike vector. 

Now, after making the definition p = n a n b T ab and recalling Einstein's 
equation, we can derive the following identity: 



Combining equations (3.1.21) and (3.1.22), we obtain the Hamiltonian 
constraint 2 : 

2 we set £ = — 1 in the last step to get the explicit form of the Hamiltonian constraint for 
a 3+1 split. Also, in the second to last line, both sides are multiplied through with a factor 



8nT ab n a n b =R ab n a n b - -Rg ab n a n b , 
8np=^(2R ab n a n b -ZR), 
16np =2R ab n a n b - £R , 



2R ab n a n b 



lQnp + ^R. 




of £, and all factors of £ 2 



are set equal to 1. 
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R=R-£ (levrp + £R) + f (K 2 - K a ^ afe ) , 
16^p = -R + ^(-K ab K ab + K 2 ) , 
16np = -£R- K ab K ab + K 2 , 

lQ-np —R — K ab K ab + K 2 . (3.1.23) 

Now, using a slightly different technique, we can derive the Momentum 
Constraint. First, project the Einstein Equation onto n a ^ b C) and define j a = 

-T 6c n c 7 b a : 

87rT ab n a 7 6 c =R ab n a -f b c - -g ab n a ^cR, 
-8tt7' c =n a -f b c g de R dbea - ^n b -f b c R, 
=2 7 V e V [rf V 6] n e , 

=2 7c fe V [d V 6] n d = 7c V,V b n d - 7c 6 V b V rf n d , 
=V d ( 7 fc c V fe n d ) - (V dl b )V b n d + V C K , 
= - W d K d + £n b (V d n c ) W b n d + £n c (Vy ) V 6 n d + V C K , 
8nj c =W d K d c + ^ cd n b V fe n d + in c K d b W b n d - V C K , 

=g bd W d K bc - in b n d W b K cd - £n c n d V b K b d - V C K , (3.1.24) 
= (g bd - in b n d ) (g c e - £n c n e ) V d K be - V C K , (3.1.25) 
=V b K b c - V C K, 
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S7ij a =V b {K ab - ^ ab K) 



(3.1.26) 



Here, Equation (3.1.26) is the momentum constraint, and in moving 
from equation (3.1.24) to equation (3.1.25), we used the fact that n a n b n c ridV a K bc = 
n a n c naV a (n b Kbc) —n a n c Kf lc ndV a n b = 0— 3 . We have now derived expressions 
for the Constraint equations for the standard ADM 3+1 split. 



3.2 Evolution Equations 

Now, we will derive expressions for the ADM evolution equations. Once again, 
rather than going with the more elegant Hamiltonian framework, we are in- 
stead going to explicitly evaluate the time derivative (in this case, the Lie 
derivative along the vector t a defined in equation (2.1.26)) of our evolution 
variables / y ab and K ab . First, let us find the evolution equation for j ab : 



£tlab =£anlab + £plab , (3.2.1) 

=an c V c 'j ab + 7cf>V a (an c ) + 7a C V 6 (an c ) + f3 c V c ^ ah 

+ 7 C 6V a /3 c + 7 a C V fe /3 c , (3.2.2) 

=a[n a n c V c n b + n b n c V c n a + V ' a n b + V b n a ] + 2V (b /3 a ) , (3.2.3) 

=a[{n a n c + g a c )V c n b + {n b n c + g b c )V c n a ] + 2V (6 /3 a) , (3.2.4) 

= - 2aK ab + 2V [a (3 b) . (3.2.5) 



3 Note that the common factor here is included simply because it also appears in the 
derivation of Equation (3.1.26) 
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Since it is expressed entirely in terms of 3-dimensional quantities (namely, 
we've eliminated any reference to n a ), we can consider equation (3.2.5) to be 
the evolution equation for 7^. Note that no reference to the Einstein Equa- 
tion was required in this derivation. We can consider this evolution equation, 
therefore, to be more of a definition of the relationship between t a , K a \, and 
a than anything else. 

Now, let us derive an expression for the evolution equation for K a b. 



£tK a b =£ an K a b + £pK ah , (3.2.6) 

=£ an K ah + (3 c V c K ab + K cb V a (3 c + K ac V b (3 c , (3.2.7) 

(£ t - £ p )K ab =an c V c K ab + K cb V a (an c ) + K ac V b (an c ) , (3.2.8) 

=a[-n c V c (Y a 7b V d n e ) + K cb V a n c + K ac V b n c ] , (3.2.9) 
=a[-n c ^ e b V c V d n e - n c : (V 'dDl^ > e 

- ™ c (V c 7a) V<m b - K cb K c a - K cb n a V c ln(a) 

- K ac K c b - K ac n b V c ln(a)] , (3.2.10) 
=a[-n c ^R cd Jn f - n c ^V d V c n e 

+ n c n b {V c n e )K ae - n c n d (V c n a )V d n b 

- n c n a {V c n d )V d n h - 2K ac K c h \ - n a K c b V c a 
-n b Ky c a, (3.2.11) 

[(7 C/ - 9 cf )lHR d cef - 7a7 6 e V,(n c V c 
+ lalb(V d n c )V c n e + n b K ae V e ln(a) 
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- (VJn(a))Wn(a) + n a K db V d ln{a) - 2K ac K c b ] 
-n a K%V c a-n b K a c V c a, (3.2.12) 

r I nk m T3 n f? 

[ /a 7 76 t^lnmk la lb -^mn 

- 7aVV d (V e Zn(a)) + ^a C ^ - 2K ac K c b ] 

--(V a)V 6 a, (3.2.13) 
a 

=a[ 7 a'7 nfc 7r^n mfc - 7a m V(8vrT mn + ^ o6 ) 

- K ac ir 6 ] - -(V a «)V fe « - aV a (-V 6 a) , (3.2.14) 

a a 

=a[R ab - K ac K c b + KK ab - 8nS ab - ^R^ a b 

-K ac K c b ]-V a V b a, (3.2.15) 

=a[R ab - 2K ac K c b + KK ab - 8irS ab + An lab (S - p)\ 

-V a V fe a, (3.2.16) 
£ t K ab =a[R ab - 2K ac K c b + KK ab - 8nS ab + 4n lab (S - p)\ 

- V a W b a + (5 c W c K ab + K cb S7 a p c + K ac V b (3 c . (3.2.17) 

In going from lines (3.2.14) to the next line, the definition T a b7 a m 7 fc ra = 
S mn was made, and the result from equation (3.1.19) was used. In going for- 
ward from (3.2.15), we contracted the Einstein Equation on g ab in the following 
way: 



g ab R ab - l -g ab g ab R =87rg ab T ab , (3.2.18) 
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R-2R =87r( 7 ab - n a n b )T c 
-R =8tt(S - p) . 



(3.2.19) 
(3.2.20) 



Since, once again, equation (3.2.17) depends only on 3-dimensional 
quantities, we can consider it to be the proper evolution equation for K a f,, and 
we now have the complete set of ADM equations, namely equations (3.2.17), 
(3.2.5), (3.1.23) and (3.1.26), and this concludes our "brute force" derivation 
of the ADM equations. 
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Chapter 4 
Hamiltonian General Relativity 



Now, we above performed the brute force derivation of the ADM equations 
for a few reasons. First, there was the goal of seeing true completeness to our 
derivation-there was no direct reference to non-geometric quantities above- 
we were able to show that the ADM equations are a direct consequence of 
merely the Einstein Equation and the rules of Lorentizan geometry. Second, 
the above procedure generalizes quite nicely to a foliation of a spacelike two 
surface embedded in a spacelike three surface-the relevant situation to many 
initial value boundary problems in General Relativity. 

This is not to say, however, that the more standard Hamiltonian ap- 
proach to this derivation is without merit. It is only through this approach 
that we can investigate several important issues-in particular, the Hamiltonian 
formalism enables us to investigate the angular and linear momentum stored 
in a spacetime, in addition to that spacetime's net energy. The brute force 
approach given above necessarily only deals directly with bulk terms-these are 
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things that are best viewed through an action principle. So, here, we will work 
through this action principle in general relativity. 

4.1 From the Lagrangian to the Hamiltonian 

We begin with the well-known (see, for instance, [10]) Hilbert action: 

Here, g a b is the four-dimensional Minkowskian metric of the spacetime, 
g is the determinant of g^, and R is the trace of its Ricci curvature. We wish 
to take this Lagrangian density and form a Hamiltonian density with proper 
equations of motion. It is a known result that no second time derivatives of 
the g ta appear in the action (4.1.1), and therefore, if we were to take all of 
the components g a b to be our canonical variables, we would then define our 
canonical momenta to be P ab = J£-, but then we would be stuck, because 
since the Hessian 

X2c X pab 

H abcd = -^4- = 4- ■ (4.1.2) 
og a bdg c d dg cd 

would have a vanishing determinant 1 , meaning that it would be impossible to 
solve for all of the P ab in terms of all of the g a ^. In particular, in transforming 

1 To see this, if we set a and c to t, then we know that the Hessian is zero, regardless of 
the value of c and d (there are no second time derivatives of the gta)- This, in turn, then 
indicates that an entire row of the Hessian Matrix vanishes, which then implies that the 
determinant of the Hessian vanishes 
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the variables of our system from the g a b and the g a b to the g a b and the P , 
there would be a noninvertible transformation operator, meaning that degrees 
of freedom would be lost in going from the g a b to the P . The Lagrangian and 
Hamiltonian formulations of this theory would therefore not be compatible. 

In order to avoid this problem, we reduce the number of variables in 
our system to the g^, which, as shown in section 2, are exactly equal to 
the 3-components of the 3-dimensional metric 7^ if we were to do a 3+1 
decomposition using t as our slicing variable 2 . Then, using Equations (2.1) 
and (2.1.25) to define a and (3$, a straightforward calculation relates g to the 
determinant 7 of 'jab. 

g = -a 2 1 . (4.1.3) 
Then, we use Equation (3.1.22) to replace R with R, giving us: 



S = Ja^(R + £K ab K ab - £K 2 + 2£R ab n a n b ) d 4 x . (4.1.4) 

Since R manifestly has no time derivatives in it, and = \{5 a c 5b d + 
3a d 3b c ){— 2a), we can see that (4.1.4) brings us closer to the canonical form 
necessary to Legendre transform the Hilbert action: every term except for the 
last either has no time derivatives, or is expressible in terms of the cjij. So, 

what of that fourth term? Let us investigate it further: 

2 of course, if we want another slicing, we can make this choice simply by coordinate 
transforming t to some other, more general function t with a timelike gradient 
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n a n b R ab =n a n b g cd R cadb , (4.1.5) 
=n a g cd (y c V a n d - V a V c n d ) , (4.1.6) 
=n a (V c V a n c - V a V c n c ) , (4.1.7) 
=V C (n a V a n c ) - (V c n a ) V a n c - V a (n a V c n c ) + (V a n a ) V c n c , 
=V a (-£V a /n(a)) + V a (n a K) - K ah K ab + . (4.1.8) 

Before we progress further, we will quickly derive two identities. The 
first simplifies the term involving the lapse function: 



V a V a ln{a) 



=V a ( 7 a6 Wn(a)) , 
a fV b V 6 a) , 



=V 

= -7 ab V a V b a+ -(V 6 a)V a 7 



ab 



7 afe (V a a)(V b a), 



a; 
1 

a 

1 



7 a W a V 6 a + V fe aV a7 ab - -(V a a)V a a 



a 



[7 aC Va (7c 6 V 6 a) - 7 ac (V a7c 6 ) V 6 a + (V 6 a) V a7 



ab 



a 



-(V a a)V a a 



= - [V 2 a + £7 ac V b a(n c V a n 6 + n b V a n c ) 
a 

-£(V 6 a)(n a V a n 6 + n b V a n a ) - -(V a a)V a a 

a 
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(V a a)V a a 



a - - £Kn b V b a + £ 2 V h aV b ln{a) + f (V 6 a)n 6 lir 



a 



1 



V 2 a. 



(4.1.9) 



a 



Which eliminates the 4-divergences of a and replaces them with 3- 
divergences. 

The second identity simplifies the term that involves the total divergence of 
n a K 3 In the following derivation double overbarred terms like 7 will refer to 
the metric intrinsic to the boundary to the 4-dimensional manifold M. 1 



3 Here, we consider the boundary to M to be divided into three parts-the initial and final 
timeslices of our evolution-namely, the endpoints of the function r, and a part topologically 
equivalent to K x 3m. We will worry about the latter below, and take the former only as 
an integral over initial and final timeslices, which we shall neglect, as the initial and final 
states of our variation will be considered fixed. 

4 Also note that following the argument by York outlined in Appendix we must make 
a distinction between the normal to boundary in the 4-manifold and its projection onto a 
3-manifold. The former will be labeled as s a = BW a R, and the latter will be labeled as 
s a = BV a R. Note that while the two are defined by the same function R, they are quite 
distinct objects-the latter is normal to the unit timelike normal n a , and the normalization 
constants are fixed, respectively, by the 4-metric g ab and the 3-metric j ab . 
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/ 
/ 
/ 



dt J d 2 Xy/q Ks a (3 a , 
dt J d 3 x^fV a (P a K) , 
d'x^-V a (f3 a K) . 



(4.1.10) 



So, therefore, this term is also equivalent to an appropriate integral over 
dm integrated over time, or equivalently, a 3-divergence integrated over the 
4-manifold. 

Now, we can combine equations (4.1.9) and (4.1.10) with their corresponding 
terms in (4.1.8) to get: 



Now, putting (4.1.11) back into the fourth term of (4.1.4) we get that the term 
involving R ab n a n b reduces to: 



n a n h R ab = V 2 a + - V a (/3 a K) - K ab K ab + K 2 . 



(4.1.11) 



a a 




) 




d 4 xa^ [-£V 2 a + V a (n a K) - K ab K ab + K 2 ] 
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1 f o 1 

d 6 xs a n a K + 



87TG 



16ttG 



- 2^7^ - ^V 2 a 



a 



(4.1.12) 



Finally, we can substitute (4.1.12) into the Hilbert action (4.1.4) to 



obtain 



S = ^ J d 4 Xy/j [a (R-K 2 + K ab K ab ) + 2V 2 a] + ^ j d 3 xs a n a K . . 

(4.1.13) 

Where we have taken £ = — 1 again to reflect the timelike slicing of our 
spacetime in the bulk. However, the explicit inclusion of the factor of £ in going 
to equation (4.1.12) shows us that there is no sign difference in the boundary 
term, no matter whether we are evaluating the boundary term on the initial 
timeslice or on the outer boundary of our spacetime. This is fortunate, as it 
avoids any contradictions on the intersection of the initial timeslice and the 
outer boundary. 

Now, we can (due to our assumption of topology K x m for our space- 
time) remove the integral over our time coordinate r from (4.1.13), and infer 
the Lagrangian density 5 : 

5 Note that the text has replaced the normal n a with a term r a . We do this in order 
to indicate that r a is explicitly unit spacelike, as we mean for the Lagrangian density to 
be defined on a single timeslice, which has a timelike outer boundary, and therefore, this 
boundary has a spacelike normal. Also note that we have converted the boundary term in 
(4.1.13) with a divergence. 
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L = ~^ a ^ R + KabR - b - K ^ + ^^c(r c K) . (4.1.14) 

lo7rG 87TG 



Note that (4.1.14) is now written explicitly in terms of 7^ and K ab) which 
is directly related to £tlab, so we have successfully separated out the "good" 
canonical variables from the Lagrange multipliers. Therefore, we can now feel 
safe defining conjugate momenta and taking a Legendre transform. Recalling 
from the ADM equation for 7^, we have £tlab = lab = —2aK a b, we therefore 
can now define the momentum conjugate to 7^: 



U ab =^, (4.1.15) 

blab 

= ^Vl(^)(Kl ab - K ab ) + J-J-^ Vc (r c r b K ab ) , (4.1.16) 
87rG 2a 4tiG djab 

= iLc^ {Kiab ~ Rab) ■ ■ (4 - Lir) 



Where we dropped the V c (r c K) term in equation (4.1.16) due to the 
fact that since V c 7 afc = 0, we therefore have, after integrating by parts 6 : 

6 Note that while it is true that |^ = S™ b n , this is not true of -^—V c K ab . In the 
latter case, it is appropriate to interchange the functional derivative with the covariant 
derivative and integrate by parts. This originates from the fact that a formal definition of 
the variational derivative is based in taking terms proportional to differentials of functions. 
The entire functional derivative should be the multiplier of the differential of the function, 
which does not work if it is inside of a derivative. 
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J- [ d 3 x^V c (r c 7 de K de ) = (4.1.18) 
= J~ l d z x^ de V c {r c K de ), (4.1.19) 
= J~ I d 3 x^ de (K de V c (r c ) + r c V c K de ) , (4.1.20) 

hab J 

= [ rf 3 ^(7 de (V c r c )^ - V C ( 7 <V)^) , (4.1.21) 

J "lab "lab 

=0. (4.1.22) 

So, now, we have our canonical variables 7 a b, their conjugate momenta 
IP 6 , and we are left with four Lagrange multipliers in a and the /3 a . Inverting 
equation (4.1.17) to obtain a expressions for K and K ab in terms of U ab gives 
us: 



K = —Tl K ab = l —{- 1 ab Tl-Tl ab ). (4.1.23) 



An obvious consequence of this is that: 



K 2 



(4.1.24) 



K ab K 



ab 



8ttG 



(7 ab n-2n ab )( 7ab n-2ii afe ), 



(j^y (sn 2 - 4n 2 + 4ir b ru) , 
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8ttG x 2 



(4ir 6 ru - it) . (4.1.25) 

'1 ' 

Giving us the useful result: 

K ah K ah — k 2 = (^) 2 (n afe n a , - ^n 2 ) . (4.1.26) 

Knowing this' we are ready to take the Legendre transform of our La- 
grangian density L in order to define a Hamiltonian density H: 

H=W h ^ ah -L. (4.1.27) 

Remembering from equation (3.2.5) that £t"iab = lab = —2ctK a b + 
2V( a /3fe), we can then explicitly perform this procedure: 



--Il ab (~2aK ab + 2V (a ft)) - {jJ^G a ^ R + KabRab ~ K " 



+-V 2 a + J—yfyV^K) 
a o7rG 



7 And note, that if we had chosen to try and write a Canonical theory based on the g^ 
and the g a \, as our variables, this program would have failed here, as we would not have been 
able to invert equation (4.1.17) in order to get equation (4.1.23), and therefore, we would 
not have been able to write down a Hamiltonian in terms of the g a b and the IT ab 
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—a 



( 



167rG 



) 




2 



) 



«\/7 
167rG 



i? + 2n a6 V (Q /3 b) 



V 2 a- 



1 



v/7V a (n a K) . 



(4.1.28) 



16ttG 



8ttG 



Which is the ADM Hamiltonian density, as found originally by Arnowitt, 
Deser and Misner [11]. A 3-dimensional integral of this quantity would then 
yield the ADM Hamiltonian. Note that it includes boundary terms, and that 
these boundary terms must be treated properly. In order to do this, one must 
be quite careful about the proper treatment of boundary charges. In particu- 
lar, it can be shown that the the above Hamiltonian can lead to inconsistencies 
if just na'ively varied. The reader is referred to the 5 for some background on 
properly treating these boundary terms. 



Now, we take the variation of the Hamiltonian density in order to find equa- 
tions of motion. Total divergences will be neglected in this section, but will 
be treated explicitly according to the perspective of Chapter 5 8 . Note that 
the above Hamiltonian must be modified by the addition of a series of surface 
charges in order for the overall variation to be zero. 

Having said all of this, it is now easy to derive the constraint equations. 

8 In this Chapter, we will use what we know about the bulk equations of motion in 
order to pick the appropriate boundary terms. Therefore, it is not explicitly necessary to 
carefully and faithfully preserve the boundary terms dictated by the Hilbert action. Instead, 
we faithfully maintain the appropriate behaviour in the bulk, and then let that behaviour 
dictate what our boundary terms should be 



4.2 Constraints 



39 



The Hamiltonian constraint is derived by taking the variation of equation 
(4.1.28) with respect to a. Neglecting the boundary terms, the result is: 



{-£c)^- KabK ^ +K2 ) ■ ^ 



If it is desired, we could have added a matter density to the original 
Lagrangian, which typically would have been a function of g ab and g a b, but 
not their derivatives. Had we done this, we would have gone through the same 
(3+1) decomposition and Legendre transform process with these degrees of 
freedom 9 . This would give a well-defined matter Hamiltonian 5^, which we 
could add to the Hamiltonian given in equation (4.1.28). Making the definition 
p = \ then makes equation (4.2.1) equal to: 



°= + ^- = V / 7 



Sa Sa Sa 



(4.2.2) 

Which is the Hamiltonian constraint, as seen in equation (3.1.23). 

Similarly, we can derive the Momentum constraint, but taking the vari- 
ation of equation (4.1.28) with respect to (3 a . As we do in Chapter 5, we do 
a little bit of algebra gymnastics, and ultimately, we discard the boundary 



9 This procedure is traced in detail for a Klein-Gordon field coupled to an external grav- 
itational field in Appendix (B) 
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term, in anticipation of dealing with it in Chapter 5. Since (3 a only appears 
in a single term above, we only worry about that single term. Furthermore, 
since U ab is a symmetric tensor, the symmetrization over the derivative of j3 a 
is omitted as irrelevant. 



_ n _ ^Einstein ^Hn 



5 
5 

'Wb 



V7 



2U ab V a p t 



n ab & ) - f3 b V a 



n 



ab 



+ 



8ttG 



V a (K ab - ~f ab K) - j l 



+ 



Wb 



(4.2.3) 



Where we used equation (4.1.17) and made the definition j a = — in 
the last step. This equation is the Momentum Constraint, as seen in equation 
(3.1.26). 



4.3 Evolution equations 

Now, we follow the standard Hamiltonian procedure to derive equations for the 
canonical field variables and n ab . 10 Namely, after deriving the appropriate 

10 For simplicity, the matter fields are set to zero in this section. We have already derived 
the correct ADM equations with matter above in Chapter 3, and we are soon to be drowned 
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Hamiltonian density and solving for the appropriate constraints, we vary the 
phase space action 



s = J d 4 x (rrV - x) 



with respect to the canonical variables. This gives: 



(4.3.i; 



5S 



d 4 x 

dr / d 3 x 



(5U ab ) 7afe + Tl ab 5 (% b ) - ^- h 5W b - ¥Lfr,. 



5U 



ab 



lab ~ 



5U ab 



ab 



ab 



^ + ± (U ab 5 lab ) - 8 lab (tl ab + 



5ti 

hab 

(4.3.2 



The first and third terms in equation (4.3.2) are simply equivalent to 



the Hamiltonian equations of motion 



5U ab 



and n 



ab _ gecond 



term is a total time derivative. We can therefore do the time integral first 
over this term, giving us J d 3 xll ab 5'~y a b over the initial and final timeslices. We, 
however, assume that our system is set up in such a way that it has a definite 
set of initial conditions, which therefore means that the value of 'jab is held 
fixed on the initial and final time slices. Therefore, any variation over these 
regions must vanish identically. We can therefore drop the second term in 
(4.3.2). Therefore, deriving the evolution equations is as simple as varying the 
Hamiltonian with respect to the 3-metric and its canonical momentum, n a6 . 
This will then give us equations for U ab and 7^. Note, however, that this does 



under a large number of terms. Adding matter back in doesn't massively increase the 
complexity of what follows, but it does enough that it is not worth our trouble here. 
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not directly give us an equation for K ao . However, an inspection of equation 
(4.1.17) shows that one can solve for K ab if one knows the values of IP 6 and 
'jab 11 - Having derived the equation for K a b in equation (3.2.17), we will not 
concern ourselves with performing this transformation here. It will just be 
noted that this transformation can be performed through a relatively direct, 
if laborious, computation. 

First, we find the equation for 'jab, neglecting terms not involving U ab 



"fab 



5 



SU ab 



a 



16ttG 

"7T 



u ab u ab -^n 2 ) +2U ab v { „ 



2a 



{16nG 



^ (u ab - ^7 ab n 



-2aK ab + 2V (a /3 6) . 



+ 2V (a /3 6) 



(4.3.4) 



Where we used equation (4.1.23) in order to convert the U ao terms into 
K ab terms. Equation (4.3.4) is clearly the same equation as (3.2.5). 

Before we derive the equation for Il ab , we derive a few properties of the 

n In accomplishing this, it is absolutely useful to use the identity derived from the defini- 
tion of the inverse metric: 

£tlab = £t (lacl^ldb) , 

= {£tlac) l^ldb + lac (£tl Cd ) Idb + lacl Cd £tldb , 
= £tlab + lacldb£tl Cd + £tlba , 
-£tlab=lacldb£tl Cd ■ (4.3.3) 
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variation of quantities with respect to d'-fai,. First, we tackle the determinant 
of the 3- metric 12 : 



— ru n/ pO-nk bst 
2! ' 

= 77 ab .. (4.3.5) 

Where the last step in (4.3.5) involved recognizing two things. First, 
that the matrix 'Jns'yki£ ank £ bst is constructed entirely from the metric and the 
Levi-Civita symbol, and thus, can depend only on these two things. Second, 
that if one were to multiply r y n slki€ ank ^ bst by 7 a b, the answer would be three 
times 7. Therefore, one can conclude the answer on the last line of (4.3.5). It 
is also not too hard to work out the variation of / y ab with respect to 7^: 



8i mn 5 

= T— 7 mr 7r S 7 Sn , 
01 ab 01 ab 



5f 



5% 



ab 



^ 1 „, 1 -,))ir / 0~lrs \ sn i „,rar„, I 7 

7rs7 + 7 7 — 7 +7 Irs 



r I IIS I ' / I c / / 1 / fro I r 

C>7afe / \°lab) \0Jab 

f Xry mn \ 

7 m(a 7'>". (4.3.6) 



12 This derivation works in any dimension, not just three-in d dimensions, the Levi-Civita 
symbol has d indices, requiring d copies of the metric tensor, and making the factorial term 
Jj. Then, the rest of this proof follows in exactly the same way as this one does. 
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And we can use equation (4.3.6) to easily change the variable of the 
variation between the metric and its inverse when it is convenient to do so-in 
particular, when taking the variation of the Ricci scalar with respect to the 
3- metric. 

Now, we embark upon taking the variation of the Hamiltonian (4.1.28) 
with respect to the 3- metric 13 : 



S\ / d 3 x 



16 nG 



a ( rr 6 n afe 
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32 nG 



V7 
16ttG 



5 lab V a V b a 



32ttG j 



^7 a ^7ab - (R a % ab + l ab 5Rj) - 2/3 c U ab 5T ab c 

167rG V '/ 



(4.3.7) 



We will deal with the rest of these terms in their due course, but let 
us first consider the last term involving the variation of the Ricci tensor in 
(4.3.7). 

13 Note that the variation of the H ab V a f3t term with respect to 7 a b depends only on the 
variation of the 3-Christoffel symbol with respect to the metric, due to the fact that the 
canonical form of the U ab is considered to be the raised one, and the canonical form of (3 a 
is the one that has the lowered index. Thus, there is no dependence on ^ a b in this term. 
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4.3.1 The variation of the Ricci tensor 

Tackling this term, we start with: 



SXmcd = -J d3x Y^G ( a ^ ab SRab) ■ (4.3.8) 

Following [11], define k = for simplicity. Our goal is then to 

compute this variation of R ab in terms of the variation of 7^. To do this, 
we first expand R ab in terms of Christoffel symbols. We then remember that 
since the difference of two connections is a tensor, therefore, the difference of 
two Christoffel symbols is also a tensor, making the variation of a Christoffel 
symbol a tensor. We can then turn partial derivatives into covariant derivatives 
relative to the central metric about which we take the variation 14 : 



« icci = - j d 3 xK a^ ab 5 (d c f ab c - d a T bc c + T ab d T dc c - T ac d T db c ) , 

= - j Vld 3 x n ai ab (d c 5T ab c - djT bc c + T dc c 5T ab d + T ab d 5T dc c 

—^ac d ST db c — T db c 5T ac d ^j , 
= - J y/i d 3 x K Q 7 fl6 {V c 5T ab c + T ca d 5T db c + T cb d 5T ad c - T dc d 5T ab c 

—V a 5T bc c — T ab d 5T dc c + T dc c 5T ab d + T ab d 5T dc c 

r~ dr-p c -p crp d\ 
ac OL db — 1 db 01 aC ) , 

= - J Vld 3 x K ai ab (V c 5T ab c - V a 5T bc c ) , 

14 also, note that variations are defined in terms of partial derivatives in function space. 
Therefore, they can be freely commuted through partial derivatives 
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= - J y/l d 3 x Ka {r h 5 d c - j da 6 b c ) V d 5T ab c . (4.3.9) 

So, we have now reduced this problem to evaluating the gradient of the 
variation of the Christoffel symbol in terms of the variation of the 3-metric. 
As with almost all variational problems, we proceed by integrating by parts 15 . 



Mmcd = ~ j Vld'x {nV d [a (r b 5 d c - J ad S b c ) 5T ab c ] 
-nV d [a{ 1 ab 5 d c - 1 ad 5 b c )]5Y ab c } , 
= - j^qd 2 xna ( 7 afe f c - f a 5 b c ) 6T ab c 

+ J ^d z xn ( 7 ab V c a - 5 b c V a a) 5T ab c , (4.3.10) 
= j^dd 2 xB ab c 5t ab c + J ^d 3 xT ab c 5T ab c . (4.3.11) 

Where we make the definitions B ab c = na ( K \f a 5 b c + \ f b 5 a c — 7 afc f c ) and 
T ab c = k (7 afe V c a - \5 b S a a - ^ a c V fe a) 16 . Note that both objects are true 
tensors whose indices can be freely raised and lowered with ^ ab and its inverse. 
Now, our variational problem has been reduced to finding the variation of a 
three index tensor multiplied by 5T ab c . In order to do this, we quickly work out 

a property relating the variation of the metric to the variation of its inverse: 

15 Note that if the factor of a were not present, we could simply pull through to 
the beginning of the equation and the integration by parts would produce nothing but a 
boundary term 

16 Since no allowed variation of the metric will induce torsion in the connexion, we know 
that (5r a {,° is symmetric on a and b. Therefore, we can symmetrize on these indices of B ab c 
and T ah r 
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h ab (7 ac W 6 ) , 

= ih ac ) i cd i db + 7 ac {hc d ) i db + i ac i cd [h db ) , 

=25 1 ab + 1 ac 1 db Si cd , 
5i ah = -i ac i dh 5 lcd . (4.3.12) 

Having this identity in hand, we now work out the general expression 
for a tensor multiplying a variation of a Christoffel symbol. Using the decom- 
position the the Christoffel symbol into metric components, and converting 
partial derivatives into covariant derivatives, we obtain the following: 



M ab c 6T ab c = l -M ah J [j cd [d albd + d blad - d dlab )] , 

= - l -M ab c ( 7 ce 7 d/ 57e/) (dalbd + d blad - d dlab ) 

+ l -M ah cl cd (dj lbd + d b 5 lad - d d 5 lab ) , 
= - M abe T J5 le} + l -M ahd (V a 5 lbd + f ab c 5 lcd + f ad c 5 lbc 

+^bhad + ^bahcd + ^bd C hac ~ V ' d &lab ~ ^dahcb - ^db^lac) , 

= l -M abc {V a 5 lbc + V b 5 lac - V c 5 lab ) , 

=- (M abc + M bac - M cab ) Vahbc ■ (4.3.13) 
Using the above definitions and symmetrizing over the indices b and c, 
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since the multiplier is V a 7f, c : 



B abc + B bac _ B cab =K(y | l^&c + -f*^ - 7 <V + -f^ + -f^ 



_^ab-c _ l^c^ab _ l^a^fec _|_ ^ca-l 

=na (f a j bc + 2f 6 7 ac - ^r c 7 afe 
(B abc + 5 6ac - 5 ca6 ) s?/mfcc =m (f a 7 fec - 7 a V ) ) , (4.3.14) 

T abc + T 6a C _ ycafe = R / yriy^ _ i^y^ - -7^^ + 7^0 

V 2 2 

2 ' 2 ' 

+ ^7 ab V c a + ]p cb V a a 

=« 07 ab V c a - 7 bc V a a - 2 7 ac V 6 a 

^ T a6c + T 6ac _ ycaft) ^ = fi . ^(fty^Q - ^^a) . (4.3. 15) 

Thus, equations (4.3.14) and (4.3.15) give us the final expressions for 
the terms multiplying V a 7h c in the variation in equation (4.3.11), since equa- 
tion (4.3.13) tells us how to convert the variation of a Christoffel symbol into 
the gradient of the variation of the 3-metric. Putting all of these equations 
together, and then integrating by parts 1 ', we obtain: 



17 note that we automatically apply Gauss's theorem to the divergence term originating 
from the bulk term. 
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+ j ^7d 3 z^( 7 a(6 V c) a- 7 6c V a a) V a 5 lbC) 

-^ 76c V a [a (fV c - 7 a( V))] + |f a ( 7 a(fc V c) « - 7 fec V a «) <5 7fec } 
+ / VT^fv a ( 7 6c V a a - 7 a(b V c) a) 5 T6c . 

= j^qd 2 x^ { V a [a (f a 7 fec - 7^) 5 7bc ] 

+5 7bc [a ( 7 bc V a f a - V ( ¥ c) ) + 2f (6 V c) a - 2 7 6c f a V a a] } 

+ J V^d 3 x^ ( 7 6c V a V a « - V (b V c) «) <5 7bc • (4.3.16) 

Now, equation (4.3.16) nearly has the form that we want. Most of the 
terms are direct multipliers of the variation of the 3-metric 7^. We are left, 
however, with a 3-divergence to be integrated over the boundary. The na'ive 
thought would be to simply eliminate this term using Gauss's theorem, since 
it should be convertible to an integral over the boundary of the boundary, and 
a well-known theorem of topology tells us that the boundary of the boundary 
is zero for any manifold. 

Finally, remembering our promise to wait until Chapter 5 to deal with 
the boundary terms, we drop everything in (4.3.16) except for the bulk term, 
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and our final answer is: 



J d 3 X1 ab 8R ab \ 7 = J d 3 x^ (Y b V 2 a - V a V b a) 8 lah . (4.3.17) 
4.3.2 The equation for tl ab 

Now, having gone through the work of that extensive aside, we are ready to 
return to the equation for Il ab that we obtained from the ADM Hamiltonian. 
Substituting (4.3.17) into (4.3.7), and using (4.3.13) to deal with the term in 
(4.3.7) involving the shift vector 18 , we get: 



5H\ 
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d x 



32nG 



87rG 



rcdi 



-IT 



32 nG 



v/7 
32nG 



a5 lab I w c w b - ^n a6 n 
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V7 

16 7TG 



8 lab V a V b a 



^ 1 (n afe /3 c - n bc /3 a - u ca /3 b ) v c 5 7 , 

87rGa (\ 



ab 



d 3 x 5^ ab 



7 afc n 2 _ 7 ^ n cd n cd + 4n a c n c6 - 2ir fc n 



4 16ttG 
+ y/j8-f ab V c 



aR ab - a-Y R ~ V a Va 



(ll ab /3 c - u ac p b - n bc /3 a ) 



(4.3.18) 



18 Once again, neglecting boundary terms, as they will be treated with in Chatper 5 
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Finally for all of that effort, we get our final solution, after referring to 
our original phase space action variation in (4.3.2): 



rr 



Slab 

V7 V 




16 



irG V 



(4.3.19) 



Equation (4.3.19) can be shown, through a laborious calculation involv- 
ing (4.1.23) and (4.3.4), to be equivalent to (3.2.17). We have now completed 
our derivation of the 3+1 Einstein equations in the bulk using the Hamiltonian 
formalism. Now, we move on and work with the boundary terms in the next 
Chapter. 
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Chapter 5 

Boundary Charges and 
Counterterms 

5.1 Gauss's Theorem 

Now, let us consider the way to extract information about boundary charges 
out of the ADM formalism. For our current purposes, we will not concern 
ourselves with the J y/jd 3 xK terms over the initial and final slices in (4.1.28). 
Instead, we focus on the term that has the form f dr f d 3 x- s / ? f'V 2 a. Before 
we progress, we will restate Gauss's theorem as described in Hawking and 
Ellis [12]: 

First, consider the integral of the following quantity over an arbitrary 
region of a n-dimensional manifold endowed with metric tensor g a b, and 
metric-compatible connection V a : 
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d n x^\g\V a V a . (5.1.1) 

Since the exterior derivative acting on the volume element necessarily 
yields zero (there are no totally antisymmetric (n+1) dimensional objects over 
a n dimensional manifold), we can apply integration by parts to the above 
term. This reduces the dimension of the volume element by one, and gives us 
an integral over the boundary of 

[ d n x^/\g~\V a V a = [ d n - l x^\r a V a . (5.1.2) 

Where the integral in (5.1.2) is over the boundary of and r a is the unit 
"outward" normal to this boundary. If one prefers to think of n dimensional 
quantities instead of n — 1 dimensional quantities, one can define the boundary 
of ^ as the zero of a scalar function p 1 . We can then define "r" = V a p. 

"r", however, differs from r a merely by a factor of ^==^. A simple 
computation then shows that (^VapV^Xv^M") = \Ad~\- Therefore, if one 
would prefer, (5.1.2) could be rewritten as: 

/ d n x^\g\V a V a = I d n - x x^\gyW" a V a . (5.1.3) 

J%. Jd%, 

For a great many applications, this conversion greatly simplifies the 
calculation of an integral, as the domain of dependence for the problem has 
been reduced from the entire bulk of the problem to merely its boundary. 

1 For most cases, we would take this to be something supremely simple, like p = r — ro, 
but it is not absolutely necessary to do so. 
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Furthermore, since the boundary of the boundary is zero by a famous theorem 
from topology, we can immediately conclude, that the right hand side of (5.1.2) 
vanishes if we are able to repeat this procedure. We should not, however, take 
this to mean that the integral of a n-dimensional Laplacian over the bulk 
always vanishes. First, we will see a simple counterexample to show that 
na'ive double application of Gauss's Theorem results in nonsense, and then it 
will be explained why this is the case. First, consider the differential form of 
Gauss's Law: 

V a E a =-p. (5.1.4) 
eo 

It is a simple matter to integrate this equation over space on both 
sides, and then apply Gauss's theorem in order to obtain the integral form of 
Gauss's Law (for clarity, we absorb the metric determinant factors into the 
volume elements): 

/ d 3 xV a E a = [ d 2 xr a E a = [d 3 x^ = ^. (5.1.5) 

J<M JdM J e e 

Now, one might then recall that the Electric field is merely -1 times the 
gradient of the electrostatic potential <f), and then replace E a with — V a in 
(5.1.5): 



Qi 



eo 



d 2 xr n V a c 



55 



= / d 2 xctN a r a - [ d 2 xV a (0r a ), (5.1.6) 

JdM JdM 

"=" [ d 2 xctN a r a - [ d 1 xs a (pr a , (5.1.7) 

JdM JddM 

"=" / d 2 x(fN a r a . (5.1.8) 

JdM 

Where we used the metric compatibility of V a to raise and lower in- 
dices both inside and outside the gradient. A computation of (5.1.8) for even 
the simplest examples from Electrostatics will show that (5.1.8) is nonsense, 
however. For instance, if we simply consider the static point charge field from 
Coloumb's law, we have = j^r^) while for r a normal to a sphere in Euclidean 
space, we have V ' a f a = - Therefore, the right hand side of (5.1.8) comes out to 
J Q 27r Jq d(j)d6r 2 sin(fl) 27re 9 or 2 , which then causes (5.1.8) to yield = ^, which 
is clearly incorrect. 

What is the origin of this error, however? It turns out that if you com- 
pute the second integral in (5.1.6), you will find that it has value which 
accounts for the discrepancy between the left and right hand sides of (5.1.8). 
But, then, what of our application of Gauss's Theorem here? It must have 
been incorrect, as it assigned a zero value to something that was not zero. 
This, in fact, was our problem, and this is the reason why care must be taken 
with Gauss's Theorem. The reason why we cannot apply Gauss's theorem to 
(5.1.6) lies in the type of derivative operator that appears in (5.1.6)— that V a 
is the derivative operator that is compatible with the bulk 3-metric, but not 
necessarily with the boundary two metric-in fact, since areas increase as we 
increase the value of r, the boundary two metric cannot be compatible with 
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V a -the bulk derivative operator should be able to record information about 
the increase of areas. Gauss's Theorem would only be applicable to (5.1.6) if 
the derivative operator appearing there were the two connection compatible 
with the two metric. Therefore, care must be taken at all times with these 
operators, lest unwanted errors creep into calculations. 



5.2 Boundary Charges in Maxwell Theory 

Now, let us extend this calculation to the full Maxwell Hamiltonian. We start 
with the Maxwell Lagrangian with arbitrary external current j a : 

L=^F ab F ab + fA a , . (5.2.1) 

here, as is usual, F a b = d a Ab — dbA a , and indices are raised and lowered using 
the Minkowski metric rj a b. Due to the antisymmetry of F a b, it should be clear 
that there are no time derivatives of A t appearing in this action. Therefore, 
we should treat A t as a Lagrange multiplier here. Now, we will 3+1 split the 
Maxwell Lagrangian along surfaces of constant t, defining A a = (0, Aj) and 
j a = {p,f). We wish to define a Hamiltonian from the given Lagrangian, using 
Ai and its canonical momentum II 1 as canonical variables. First, we solve for 
IP: 



n ' s H = ^rni-f-" = F " = d ' A ' - d ' A ' ■ < 5 - 2 - 2 > 
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So, now, the Maxwell Lagrangian can be rewritten as: 



L= l -F ah F ab + fA a , 

= l -F u F u + l -F^F l3 + fA t + j% , 
= - + \f'^F 13 + p0 + fA t . (5.2.3) 

Now, we perform the Legendre transformation to arrive at the Maxwell 
Hamiltonian density: 



X =11% - L , 

1 • 1 

2 1 4 



rr (-n, + dtf) - ( --u% + -F^F tJ + p( p + fA t ) , 



- - \f^F 13 + 11^0 - P 4> - fAi . (5.2.4) 



Now, let us take a variation of the phase space action with respect to 

0. 



5\ 4> X = J d 4 xirdi5<i)- p5(j>, 

dt J d 3 xdi (IT0) - S^diW - P 5(p , 



(5.2.5) 
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the term under the J d 3 x is, upon realization that IP = — E l , easily recog- 
nizable as Gauss's Law. The boundary term, however, is the integral form 
of Gauss's Law, only without an appropriate charge term to balance it. Fur- 
thermore, the value of cf) in this formalism is pure gauge-there is no abstract 
reason it has to take any value at all. Therefore, the presence of the boundary 
term gives us no reason at all to believe that the action derived from (5.2.4) 
is stationary with respect to variations of 0. 

Therefore, in this theory, the behaviour of the bulk can create problems 
with the maintenance of a stationary action, by the production of gauge- 
dependent boundary terms. This behaviour, however, can be controlled. 

Consider an ansatz to resolve this inconsistency-we define the Hamilto- 
nian as the 4-dimensional integral of the Hamiltonian density plus a boundary 
term dependent only on and an as-yet undetermined constant: 



1 



1 ™-„- 



H = d x — IHL - -F^Fu + IT «9,, 



3% 



+ d> d' z xQ(j)+ (t d 2 xJ i A i 



(5.2.6) 



The variation of this Hamiltonian with respect to Aj and II* can be 
shown, relatively easily, to give Faraday's law and the Ampere-Maxwell law 2 . 

2 Since there is a gauge invariance built into fy, a similar boundary term to the one 
described above is generated by the variation with respect to A;. It is dealt with in the 
same way as the variation with respect to <fi, by the addition of a set of constants J 1 that are 
determined in order to cancel the boundary integrals. Since this treatment isn't difficult, 
and is near-identical in form to the boundary terms arising from the variation of 0, these 
terms are not explicitly treated here 
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Vj-B* = is given by the fact that Fij is the exterior derivative of A4 and the 
second exterior derivative of any quantity is zero. We therefore expect the 
variation of the Hamiltonian with respect to the only remaining variable, 0, 
to give Gauss's law: 



SH _5 

50 ~50 
_ 5 
~50 
5 



d 3 x (ir<9i0 - pip) + Q j d 2 x (0) 
d 3 x (di (nV) - WT - fxp) + j d 2 x (Q(f>) 

d 3 x (djT + p) + j d 2 x (Q + nlT) , 
J d 3 x (p - ViE 1 ) + j d 2 x (Q - nBf) . 



(5.2.7) 



Since F u = IP = —E l , we recognize the integrand in the first term 
above as simply being the electric field for the given charge configuration. 
Since this equation has no time derivatives, we can therefore consider this 
first term as a constraint upon the equations of motion given by varying the 
Hamiltonian with respect to Aj and IP. Once we have solved the appropriate 
equations of motion, we can then evaluate the second term in the boundary 
integral in equation (5.2.7), and then assign the value to Q that cancels that 
boundary term. This process is not unlike the renormalization process used 
in quantum field theory-where certain terms are computed using prescribed 
Hamiltonian dynamics, and then appropriate counterterms are assigned the 
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appropriate values at the end of the calculation in order to yield consistent 
equations of motion. 

This aside being taken, we now discuss the topic of boundary charges 
in the ADM formulation. 

5.2.1 ADM Formulation 

Now, consider the ADM Hamiltonian: 



The V a {n a K) term is acted on with a 4-dimensional divergence opera- 



and final time slices, but nowhere else, provided that the spacelike normal to 
the outer boundary is normal to the timelike normal 3 . We therefore do not 
concern ourselves with this term. Now, let us do some massaging to the second 
to last term, involving the shift vector j3 a : 



3 Note that Appendix D shows that the inner product between the radial vector and the 
timelike normal to the 3+1 slices is given by a boost parameter at infinity. We can therefore 
appropriately choose a slicing that gives a zero boost parameter between the two at infinity. 
In the asymptotically flat case, this is equivalent to the condition that j3 a - 




^f{V 2 a + aV a {n a K)) + 2U ab V a (3 b ] + H 



(5.2.8) 



tor. Therefore, it contributed a quantity of action J y/jd 3 xK 



over the initial 
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We group the second term above along with the first and second terms 
in (5.2.8) and the matter Hamiltonian into the generic term H&, and refer 
to them as the bulk Hamiltonian density. Substituting the above result into 
(5.2.8) gives us: 

H = Hfe + 8^G j ^v^V a (-^AJF 6 ) - J d 3 x^V 2 a . (5.2.9) 

Now, since we have two integrals of divergences over the bulk manifold, 
we can convert these integrals into integrals over the boundary of the spacelike 
slice: 



H = Hb+ 8^G J d2x V^ r Ml ab K-K ab )-^J-^ J d 2 x^r a V a a. (5.2.10) 

Now, we might be tempted to deal with the V a o; term by pulling the 
V a through, integrating by parts, reapplying Gauss's Theorem, and relying 
on the vanishing of the boundary of the boundary. Hopefully, the example 
from Equation (5.1.2) provides sufficient intuition to show that, since we have 
a V and not a V appearing in (5.2.10), this path will lead us inexorably to 
nonsense. Instead, we must consider the Boundary Hamiltonian to have the 
form provided above, and to take this seriously as we take our variations, which 
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we will do below. 

5.2.2 Variation of the Boundary Hamiltonian 

When we take the variation of the Hamiltonian with respect to the lapse 
function, we obtain: 



9H dH b 1 f 9V a a 



da da lQnG J ' da 

1 r* 



=C ~TQ^Gj d xV ^r a )- (5-2.11) 
Taking ^ simply gives the Hamiltonian constraint 

J d 3 x^ (p - ^q( R ~ KabK -b + K 2 ) 

where p = M ™ a J ter , Now, in order to simplify this term, we are going to 
have to work through some somewhat intricate mathematics. First, we take 
the boundary to be the zero of some function /, as above. Then, we define 

P± = ^/^ /v a / ' ^' ien ' r ' a = fi-L^af' * s manifestly normal to the boundary and 
of unit length. Then, assisted by the identity ^7 = f3±y/q, we can work out: 



<9H 

da 



C = ~Y^gJ d2x ^Mr a ) (5.2.12) 
= - Yttc J d2x ^a(VvP±V a f) , (5.2.13) 
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16ttG 
1 

16ttG 
1 

IQttG 



d x 



j d 2 xV a (^V a f), (5.2.14) 
l — j d 2 x (V a /V av ^+ v^V a V a /) , (5.2.15) 
7^ a V aV /^ + ^jj ab {d a d b f - r a6 c Vj) , 

(5.2.16) 

" ma S d ' x {jS lr ° fJ + '^OAf ~ T »'k) ' 

(5.2.17) 

16^G / Sx ^~ q (~ rCf ca<1 + ^ r ^^ ~ ^ abd adbf) , (5-2.18) 



d 2 Xy/q ( ^ 7 a V c [{db'yac + <9 a 7 cb - <9 c 7 a6 ) - (<9 b 7 ac + <9 c 7 ab 



16ttG ./ v * V 2 

- <9 a 7 fe , 

1 



16ttG 



( 7 a V c [d alcb - d clab }) . . (5.2.19) 



We discarded the term ^ ah d a d b f in (5.2.18) by requiring that we define 
the surface on the boundary by a function for which this term vanishes there. 
This can always be done, since we only require that f vanish on the boundary in 
order to define /. Therefore, we can always multiply / by another function, and 
then solve the appropriate Lapace problem in order to find a new / which has a 
zero on the boundary, but which also satisfies ^ d a d b f = 0. This then shows us 
that the variation of the Hamiltonian with respect to a leaves us with only the 
term (5.2.19), which is known as the ADM Mass. For asymptotically Cartesian 
coordinate systems, it is relatively easy to show that this term will very easily 
yield the mass parameter for simple examples such as the Schwarzschild, Kerr, 
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and Vaidya metrics. 

However, it is very quickly evident that this term is in no way coordinate 
independent. In particular, if one were to simply write down the metric for 
Euclidean space in spherical coordinates, and na'ively compute the "mass", one 
would find the following, starting by computing the 3-metric: 



lab 



( 



\ 



1 

Or 2 
r 2 sin 2 6 



\ 



J 



We can define / = r — tq, which makes our surface one of constant r. 
This then yields = 1, r a = (1,0,0), and yjq = r 2 sin8. Putting all of this 
data into (5.2.19), we get the result: 



1 



Sx^q ( 7 a V [d alcb - d clab ]) = — — / d 2 x r 2 sin6 



16nG 1 UIW C,UUJ/ 16ttG 



G ' 
(5.2.20) 



To say the least, this is a silly expression for the amount of mass or 
energy contained within a sphere of radius r in the Euclidean space! In par- 
ticular, as we take the limit r — » oo, it should be clear that this expression 
diverges. Euclidean space, however, is empty, so we would expect that the 
value we derive would be zero. Furthermore, if we were to write the Euclidean 
metric in Cartesian coordinates, it should be manifest that (5.2.19) is equal to 
zero. We therefore see that this problem is, in fact, a problem with our coordi- 
nate system. Older references such as [10] and [13], and [11] will deal with this 
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apparent contradiction by demanding that formula (5.2.19) is only applicable 
for asymptotically flat metrics in which the coordinate system asymptotically 
approaches the Minkowski metric at an appropriate pace as r — > oo. Following 
work by Hawking [14], newer references (see [15]) account for this oddity by 
including a counterterm in the action equal to +^ in the case above, or, more 
generally, equal to negative one times the ADM Mass of an empty space in the 
appropriate coordinate system. This fixes this problem by making the formula 
for the ADM mass coordinate invariant, but it adds the problem of requiring 
that a fiducial flat metric be introduced into the formalism in order to define 
the counterterm. Below, an alternative technique to this counterterm will be 
introduced. It will still retain the problem of requiring a fiducial flat metric, 
but it will have the advantage of being more compactly defined, and will also 
have the advantage of being manifestly coordinate invariant. 

So, begin by defining the fiducially flat metric -y^. It should be defined 
in such a way that its coordinates are adapted to the spacetime in question-at 
the minimum, it should have the same coordinate transformation tensor as 
the true metric does. For example, if you are dealing with a Kerr solution 
to Einstein's equation, the fiducially flat metric can be defined by setting the 
mass parameter of the solution equal to zero. Then, the idea is to treat the 
true metric as if it were a tensor field living in the flat metric space. Then, 
we can define a torsion-free derivative operator T> a that is compatible with the 
fiducial flat metric (i.e., so that £> a 7b c = 0). Once this is done, we use the 
same function / to define the surface at the boundary of our space, only now, 
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we use 7 a fe to define a /3± and r a . We then do the integral in (5.2.19) over the 
fiducially flat space, replacing all of the d a terms with T> a terms. The final 
result is: 



M ADM ,ne W = J V^d 2 x^f ab r c {V albc - <D clah ) . (5.2.21) 

This expression has several advantages. In the case of asymptotically 
flat coordinate systems, all of the r a (, c 's are all zero, the metric is equal to 
^ab + O(-) at the surface at infinity, and therefore, it should be clear that 
(5.2.21) gives the same result as (5.2.19). It should also be clear, moreover, that 
the integrand of expression (5.2.21) transforms as a scalar under coordinate 
changes that affect the fiducially flat metric, provided, as stipulated when we 
defined j a b above, that both metrics have the same coordinate transformation 
tensor. Therefore, it should be clear that (5.2.21) makes sense if the bulk 
metric is defined in any coordinate system, so long as that coordinate system 
covers an appropriate two-surface at the spacetime's spacelike boundary. 

Furthermore, as a concrete example, consider the Kerr spacetime in 
the spheroidal coordinates given in Appendix A. Here, the 3-metric is directly 
obtainable from the 4-metric, and we can find the fiducial metric simply by 
setting M = 4 : 

4 As always in this work, defining A = r 2 + a 2 and B = r 2 + a 2 cos 2 9, and using the 
spheroidal coordinates defined by taking x — yAsindcosip, y — \fA~sin6 simp . and z = 
rcosip, followed by a second coordinate transformation d(f> = dip — \dr. Coordinates are 
then rendered in the form (r, cp, 9) . This form for the metric tremendously simplifies the 
expression for the metric when compared to the original form, which is a dense 4x4 matrix 



67 



/ 



lab 



2Mt- 



"(1 + 



2Mr- 
B ■ 



asin 2 9 



;i + ^L)asin 2 6 



B 



^) sin 2 9 



lab 



V 



o 

1 —asin 2 6 

—asin 2 6 Asin 2 6 

£ 



5 



\ 



/ 



(5.2.22) 



Then, it is simple enough to compute the inverses of these metrics, and 
follow through with the standard prescription given above, and to compute 
(5. 2. 21). 5 Doing this gives us f a = (y^f, 0,0), and then it only requires the 
basic definitions of the covariant derivative to compute (dropping the factor 
of G now, since the metric is written in geometricized units): 



,new ) 

1 f 2 f MsinO [16r 4 + 4a 2 r 2 (cos(2fl) + 3) + a 4 (cos(4fl) - 
~\M J X \ 4 (r 2 + a 2 cos 2 ef J 

(5.2.23) 

The integration, of course, is easy. All that it does is cancel a factor 

of 27r in the numerator against the factor of 16ir in the denominator. The 9 

in both its covariant and contravariant forms. 

5 Of course, we are free to compute (5.2.19), too. Since, however, this metric reduces to 
a spherical expression of the Schwarzschild Metric for a = 0, we would find that this would 
be beset with the exact same singularity discussed above 
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integration can be computed with a package such as Mathematica, yielding, 
amazingly, simply M, independently of which radius is chosen for the inte- 
gration. This indicates that (5.2.21) has a flavour much more similar to that 
of Gauss's Law than the expression that you would get from (5.2.19), which 
would only give you the correct value of M if you were to do the integral on the 
sphere infinity with the metric expressed in the original Kerr coordinates 6 , or 
another asymptotically Cartesian form. This technique, however, is both man- 
ifestly coordinate invariant and works at an arbitrary radius 7 without adding 
explicit counterterms to the action. Furthermore, it automatically yields a 
zero result for an empty space, simply by the construction of <D a . The only 
disadvantage is that it required a massaging of terms after taking a variation, 
rather than explicitly coming out of a variational principle, in the manner that 
the terms derived in [11] and [15] were. It also retains the problem, inherent 
to all of these approaches, of defining some sort of flat metric with which to 
compare the true metric, whether it is done implicitly by choice of coordinates 
in the style of the original ADM Paper, or more explicitly, as shown here or 
in the counterterm approach. 

We can apply a similar procedure to the shift vector /3 a . Taking only the pure 

gravity terms in the ADM Hamiltonian involving /3 ffl , we have: 
6 see footnote number (4). 

7 while this slicing is invariant under a 3-dimensional coordinate transformation, and al- 
ways yields the ADM mass at conformal infinity, at other radii, it is not slicing independent- 
this alternate ADM mass, when computed in a Boyer-Lindquist slicing, gives simply M at 
infinity, but will monotonically increase to infinity as one approaches the horizon 
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H involving p = (- f d*x^p b V a ^W b + j d 2 x^^r a f3 b W b ^j . 

(5.2.24) 

Varying this portion of the action with respect to /3 a , and remembering 
that IT 6 = ^7 (K ab - j°*K) gives us: 



6H 



5(3 a 8ttG 



J d?x^V a (K ab - ^K) - j> d 2 xj-qr a (K ab - j ab K) 

(5.2.25) 



The first term above is easily recognized as the 3-dimensional integral 
of the Momentum constraint, which must therefore vanish. The second term 
above, however, does not, in general vanish. In particular, it can be shown to 
be equal to Mv % for boosted black holes, and equal to M aip 1 for Kerr black 
holes. Therefore, this term cannot be expected to vanish for a general solution 
to Einstein's equations. Just as in the Maxwell case above, and for the charged 
case, however, this can be treated simply by adding a boundary counterterm 
§ d 2 Xy/q (—P a /3 a ) to the Hamiltonian. Then, after the bulk equations of mo- 
tion have been solved, we can then solve for the appropriate value of P a such 
that it cancels the second term in the variation given in equation (5.2.25), 
thereby giving an overall variation of the Hamiltonian that is equal to zero. 
Note that each of these counterterms are, in fact, charges corresponding to 
conserved quantities arising from symmetries in the Hamiltonian/Lagrangian. 
The Lapse function is associated with time reparameterization invariance of 
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General Relativity, and therefore, corresponds to an asymptotic energy. Mean- 
while, the shift is associated with the choice of 3-dimensional coordinates on 
the 3+1 fibers. Therefore, the shift vector encodes information about asymp- 
totic translations and rotations. Therefore, the charge, P l associated with it 
can be associated with linear and angular momentum. 
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Chapter 6 
Null Geometry 



In the following chapters, we are going to be investigating several properties 
of null submanifolds of general spacetimes. Null spaces have quite a few spe- 
cial properties that are somewhat nonintuitive if one is used to dealing with 
nondegenerate vector spaces. In particular, there is no natural connection be- 
tween the intrinsic vector space tangent to a null manifold and the equivalent 
covector space. Consequently, one must take care in constructing the induced 
geometry from an enveloping spacetime. For this reason, explicit examples are 
worked out in Appendix A 

First, remember the techniques derived for projecting onto a surface 
derived from section 2, summarized here for convenience: 

1. Find a function / such that / = constant on the surface in question 

2. choose / as a coordinate in the spacetime in question 

3. Eliminate the df terms from all relevant forms, and in particular, the 
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lowered version of the metric tensor. Call this new metric tensor the 
induced metric on the surface 

4. Take the inverse of the induced metric, and use this as the raising oper- 
ator on the surface 

5. If there are vectors that need to be projected, you can use one of the 
following (equivalent) techniques: 

(a) If the vector already satisfies v a df a = 0, then the appropriate zero 
component can simply be dropped, and the vector can just be con- 
sidered to be a vector living in the fiber. Otherwise: 

(b) lower the vector using the 4-metric, use the above procedure to 
project, and then raise using the inverse three-metric 

(c) define n n = , Va ^ =. Raise n„ using the four-metric. Then 
project vector indices using the operator ^ a b = 5 a b —!;n a n b , as v ar ) a hl 
will, by construction, have zero components along the appropriate 
direction, which can then simply be dropped. 

6. If it is necessary to project back into the 4-dimensional space, 3-dimensional 
vectors can simply be included with the appropriate zeros added 

7. If you wish to project one-forms back into the 4-dimensional space, how- 
ever, it is necessary to operate on them using the 7 a 6 operator 

1 Recall that £ is the sign of n a n a 
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Now, we are about to need to do the above for the case of a null subman- 
ifold of a 4-dimensional space. It might seem trivial to just follow the above 
steps that supposedly would allow us to seamlessly go back and forth from 
the 4-space to the 3-space. But on closer glance, step (iv) requires that we be 
able to invert the induced 3-metric. If the 3-space is null, then the 3-metric 
necessarily will have an eigenvector with zero eigenvalue, and therefore, will 
not be invertible. Also, in the null case, the vector n a defined in (v)(b) will 
be singular, as V a /V a / = 0, which means that the projection/inclusion oper- 
ator 7 a 6 cannot be defined. Therefore, in the null case, some extra care must 
be taken in the above procedures when working out the appropriate induced 
geometry. 

So, first, we start by noting that, until the fourth step, there was nothing 
problematic at all. In particular, we still consider our surface to be a solution 
set for some equation. In particular, for a Kerr spacetime in spheroidal Kerr 
coordinates, we consider the horizon to be the solution to the equation r a + 
a 2 — 2Mr = 0. This makes the first few steps easy enough to follow. In 
particular, this means that we can map forms in the 4-dimensional spacetime 
to forms on the null sub manifold using the same procedure (i.e., by dropping 
the r component of the form in the Kerr spacetime). Note, however, that 
this technique will already create some anomalies for null sub manifolds. In 
particular, it causes the lowered version of the outgoing null normal to the 
Kerr horizon to be projected to the zero vector on the horizon. This might 
lead to some mystery regarding what the three covectors that span the dual 



74 



tangent space to the Kerr horizon are, but before this mystery is settled, let 
us soldier along. 

First, let us denote the 3- metric by q a b 2 ■ If the form that we omitted in 
order to define q a ^ was null according to the 4-metric, then it is the case that 
q a b can be shown to have zero determinant. This then means that there is a 
three vector that is mapped to the zero form by q^. To ascertain the nature 
of this vector, note that, by the fact that df a is null, that df a g ab will have no 
component along the df direction. We therefore can infer that the vector that 
is mapped to the zero covector by q a b is, in fact, g ab dfb, times a function, which 
we will henceforward refer to as the outgoing null normal l a . 

Now, let us return to looking at the basis of the covector space. If 
we use q a f, to evaluate the norms of the vectors in the tangent plane to the 
null surface, we can see, relatively quickly, that the tangent plane to the null 
surface is spanned by i a and two other spacelike vectors. Therefore, we would 
expect to have the covector space spanned by three vectors that, when acting 
upon unit vectors in the tangent space, will give you Kroneker delta functions 
as output. Since the raised versions are already spacelike under the action of 
q a b, for the spacelike vectors, we can simply raise and lower indices using the 
metric of the horizon just like any other vector. q a ^ h = 0, however, so we 

cannot do this. And we already established that the 4-dimensional covector 

2 There is some trickiness here-it is typical practice to also define the spacelike metric 
of the 2-dimensional sections of the horizon defined by constant null parameter as q a b- We 
follow normal convention here, which is at least aided by the fact that q a b is degenerate 
with a zero eigenvalue being precisely the difference between the 2-space and the 3-space. 
Regardless, the reader is advised to be careful when interpreting something written down 
as q ab . 
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£ a has a vanishing pullback onto the null surface. However, if we specify two 
spacelike vectors and £ a , then there is a unique null covector k a such that it 
is orthogonal to the two spacelike vectors and g ab £ a h — — 1. The pullback of 

this vector onto the horizon (denoted k a ) will be the third one-form spanning 

<- 

the covector space. Then, we will require that q ab k a = 0, which guarantees 
that q ab is degenerate, as well as uniquely specifying q ab . Finally, we conclude 
this section by pointing out that we now can define the projection operator of 
vectors onto the null surface as: P a b = —k a £ b + q ac q c b, which manifestly has no 
effect upon the two spacelike vectors normal to k a and £ a , and has no effect 
on £ a , but maps k a to zero. 

Now, with both q ab and q ab uniquely specified, we can define the in- 
trinsic connection V a to the null surface according to (for v b tangent to the 
null surface) V ' a v b = d a v b + f a6 c , where f a6 c = \q cm (q am ,b + qbm,a - q a b,m)- 
This concludes our section on the embeddings of null submanifolds. We now 
progress to analyze Ashtekar's Isolated and Dynamical Horizon formalism. For 
the sake of clarity, in the appendix, we work out this decomposition for two 
null surfaces: a null cone in Minkowski spacetime and the event horizon of a 
Kerr spacetime. 

6.1 Null Expansions 

Now that we have defined our two null normals and come up with an unam- 
biguous way to define induced geometries on null surfaces, let us consider a 
few properties of these null surfaces. One quantity that will repeatedly be of 
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critical importance in the following proofs will be the notion of the expansion 
of a null normal. As was shown above, the covector space to the horizon will 
contain one null one-form, labeled k a , while the vector space will contain a vec- 
tor £ a that is not metrically related to k a in the enveloping spacetime. We will 
denote k a as the ingoing null normal, and £ a as the outgoing null normal. Note 
that £ a is the null generator of the horizon. Now, consider the null expansions 
of these two vectors, defined by 9 = q ab V ' a £ b and 9^) = q ab V a k b . These will 
have a few interpretations in this work, but for now, we will interpret them as 
the time rate of change of two surfaces along flowlines of the relevant vector. 
If one doubts this interpretation, it is easy to show that this is in fact valid: 

£ e A=£ e J ' y/q = J 2~^ £eq ' 

= \J Vqq ah £eq ab = \ J y/q^ (td c q ab + 2q ac d b t) , 
Vq (q ab t C d c q ab + 2q ab q ac (V b £ c - T bd c £ d )) , 
^ (29 {e) + q ab £ c d c q ab - q ab q ac £ d g cm (g bm , d + g dm ,b ~ g bd ,m)) , 
= l - J y/q (29 {e) + q ab £ c d c q ab - q mh £ d (g bm , d + g dm , b - g bd , m )) , 
= \f (2^(0 + q ab td c q ab ~ q mb £ d d d (q bm - £ b k m - £ m k b )) , 

= / VqM ■ (6-1.1) 

Therefore, if 9^) is constant on the two-surface, which will typically be the case 
in this work, then the rate of change of that surface's area as one Lie drags 
along £ a is simply 0(e) A, and if 0(e) = 0, then the surfaces have a constant area. 
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In the dynamical case, this argument will break down, but the reason for that 
is that £ a will cease to be a tangent vector of the dynamical horizon, which 
will be spacelike, rather than null. 

6.2 Raychaudhuri's equation 

In this section, we will derive Raychaudhuri's equation for the case of a null 
vector. This result will be used frequently throughout this text, so it is useful 
to show its origin here. 

Consider a non-affinely parameterized null vector £ a . Define its affine 
parameter by £ a V a £ 6 = n£ h . Furthermore, define a second associated null 
vector k a such that g a bk a £ b = —1, and raise and lower indices on £ a and 
k a freely using the metric tensor and its inverse. Then, the operator q a b = 
gab + £ a kb + k a £b annhilates both £ a and k a , as can be verified by inspection. 
We wish to consider q a t to, in a sense, be the metric of a two surface, but note 
that it is degenerate 3 , and therefore, has no natural inverse. We therefore must 
raise and lower indices on q a b using the full metric tensor, as is customarily 
the case when dealing with objects in a null geometry. 

Now, consider the tangent gradients of £ a given by q a c Qb d ^c^d- Obvi- 
ously, this is a tensor with two indices, and therefore, can be decomposed 
into an antisymmetric part, a traceless symmetric part, and a trace, as shown 
below: 

3 Having, after all, two linearly independent null vectors! 
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q a c qb d V c £ d = -q ab + a ab + u ab . (6.2.1) 

where the factor of | is included so that q ab V a £ b — 6- It will be of great 
interest to evaluate the evolution of quantities along flowlines of £ a , which will 
usually be taken to be the null tangent vector of the horizon. Furthermore, 
the quantity 9 will prove to be the critical central element in the definition of 
Isolated and Dynamical horizons. Therefore, it is natural to want to calculate 
the quantity £$. So, let us now do so 4 . Note that similar derivations can be 
found in several sources, including [13] [15] 

£ t e =£ e {q ab V a e b ) = q ab £eV a £ b + (V„4) £iq ab , 

=q ab £ c V c VJ b + q ab (VJ b ) V a t + q ab (V a Q V b t 

+ (V a 4) {£ C V c q ab - q cb VJ a - q ac VJ b ) , 
=q ab £ c R cab % + q ab £ c V a V c £ b + q ab (V c £ b ) V a £ c + q ab (VJ C ) V b £ c 

+ (V a 4) {£ c Vc {£ a k b + £ b k a ) - q cb V c £ a - q ac V c £ b ) , 
= - R ab t£ b + q ab V a (£ c V c £ b ) - q ab (V a £ c ) V c £ b + q ab (V a £ c ) V c £ b 

+ q ab (V b £ c ) V a £c + (V a £ b ) £ a £ c V c k b + (VJ b ) k b £ c V c £ a 

+ (V a 4) £ c k a V c £ b - q cb (V a 4) V c r - q ac (V a 4) V c £ b , 

= - R ab t£ b + q ab V a (K£ b ) + q ab (V b £ c ) V a 4 + K£ b £ c V c k b 

4 Note that g ah R ac bd c t d = q ah Racbd c t d = Rad a l h , since all of the terms by which g ab and 
q ab differ include factors of £ a , and the antisymmetry of the Riemann tensor will guarantee 
that these terms vanish. 
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v a e b + K £ b k a vj b - q ab (v c 4) v a r - q ab ( vx) v b t , 




(6.2.2) 



where the Einstein equation was used in the last line in order to convert the 
Ricci tensor into a stress energy tensor. 

Equation (6.2.2) is the null Raychaudhuri equation. It is of critical 
importance in several applications involving various horizons in relativity. In 
particular, its use was critical in the singularity proofs of Hawking and Pen- 
rose that showed that the existence of a closed trapped surface 5 along with a 
reasonable restriction on the matter content of spacetime, you will necessarily 
have a spacetime singularity. 



5 i.e., one for which, given a closed 2-surface and its two associated null vectors, £ a and 
k a , you have 9^ < and 9^) < 0. 
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Chapter 7 
Isolated Horizons 



Now, let us consider the Isolated and Dynamical Horizon 1 framework as de- 
veloped by Ashtekar and others. Previous attempts to define black holes have 
been fraught with difficulty The most typical definition of the boundary of 
the black hole was the event horizon, which is defined as the past develop- 
ment of timelike future infinity. While this definition has a convenient feature 
of capturing the nature of 'the point of no return' that is typically associ- 
ated with black holes, it also faces the significant problem that it is defined 
only globally in a spacetime. Therefore, locating a spacetime's event horizon 
requires knowledge of the entire future development of that spacetime. In 
particular, there are simple analytical examples (i.e., the Vaidya metric where 
M(v) is nonconstant for a compact subdomain of v 2 ) where event horizons ex- 

1 Henceforward, we will use the abbreviation IH to denote Isolated Horizon, and DH to 
denote Dynamical Horizon. 

2 As a concrete example, take M(v) = for v < 0, M(v) = Mqv for < v < 1 and 
M{v) = Mq for v > 1, as shown in figure 7.1. There will be regions in the past domain of 
dependence of the late-time singularity (More succinctly, points inside of the event horizon) 
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ist in flat regions of spacetime. Therefore, local observers will have no means 
by which to identify whether or not an event horizon is present. Further- 
more, beyond this conceptual 'epistemological' objection, there is the much 
more practical problem that the global definition of the event horizon makes 
it difficult to locate its intersection with a particular Cauchy slice of a space- 
time, and therefore, very difficult to locate numerically (and in fact, for the 
Vaidya example discussed above, the event horizon cannot be located using 
only quantities defined locally on that particular Cauchy slice, as the points 
for which v < were exactly isomorphic to finitely sized regions of Minkowski 
spacetime.). However, having at our disposal all of the machinery of isomet- 
ric embeddings and manifold theory in the DH case, while the behaviour of 
the horizon can still be, in the abstract, wild and superluminal, we now know 
something of how to track the motion of the DH across the Cauchy slice. 

Therefore, the focus began to be shifted toward finding a definition of 
the 'point of no return' that could be defined locally, or at least quasi-locally. 
The first attempt to do so was done by Hawking, with his notion of Apparent 
Horizon, which he defined as the boundary of a 'trapped region'. This was an 
improvement over the event horizon as a notion, as it is definable at a moment 
of time. It still proved problematic, however, as it does require some use of 
the non-local notion of a 'trapped region'. In particular, this problem would 
lead to several cases where the behaviour of the apparent horizon proved to 
be 'wild'-apparent horizons would move and evolve discontinuously in certain 

that will have completely flat intrinsic and extrinsic geometries. Observers at these times 
will have no way of knowing that they are inside of an event horizon. 
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Figure 7.1: A Penrose- Carter Diagram of a Vaidya spacetime with a mass 
function that is linear for a finite advanced time. The dashed line represents 
the event horizon, and the solid line in the interior of the diagram represents 
the dynamical horizon. 7\ represents a spacelike slice which intersects the 
event horizon, but has zero extrinsic curvature and a Euclidean metric. T2 
represents a spacelike slice that intersects the dynamical horizon. All slices 
after the dynamical horizon merges with the event horizon are identical to 
slices of Schwarzschild. 
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spacetimes for certain slicings. Still, they have proven useful, and have been 
implemented numerically. 

All of these problems, however, have been removed by the closely re- 
lated notions of Hayward's trapping horizon, and Ashtekar et al.'s notion of 
Isolated and Dynamical Horizons. These notions require merely that the no- 
tion of non-expansion be confined to a single null or spacelike submanifold 
of the full spacetime. In particular, this guarantees that the wild behaviour 
discussed above can be avoided-since the horizon is defined to be a smooth 
manifold, the behaviour of its slicings is guaranteed to also be smooth. It is 
therefore considered to be quite compatible with a great number of approaches 
in Relativity, including numerical ones. Using the machinery developed in the 
previous chapter, we now work to develop this framework in detail. We start 
with the much simpler and restrictive notion of Isolated Horizon, and move 
on to the Dynamical horizon. Finally, this chapter will conclude with a proof 
of the Laws of black hole dynamics using this formalism, and then with a 
concrete application to the Vaiyda spacetime. 

Now, consider a 3-manifold that is a subspace of a Lorentzian 4-manifold. 
Following Ashtekar and Krishnan, this three manifold will be considered to be 
a non-expanding horizon if the following three conditions are met: 

1. The 3-manifold is null, topologically R x § 2 or R x T 2 and has the null 

vector £ a as a tangent vector, and associated (degenerate) metric q a b 

2. On the three-mainifold, the expansion 9(g) = q ab V a h vanishes 

3. The 4-metric g ab satisfies Einstein's equation, and all relevant matter 
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equations of motion are satisfied, and any matter involved obeys the 
Dominant Energy Condition. 

The third condition just specifies that we are in fact talking about Gen- 
eral Relativity, and the Dominant Energy Condition guarantees the vanishing 
of certain terms in the Raychaudhuri equation, as will be seen below. Mean- 
while, the first condition guarantees that you have a static, non-evolving hori- 
zon, as translation along the null vector £ a leaves one on the three-manifold. 
The essential condition, therefore, is the second one, which captures the essence 
of the Horizon being the outermost 'point of no return': The null geodesies 
that have tangent vector £ a simply transport null observers along a horizon of 
constant area 3 , rather than allow them to move spatially, as would typically 
be the case for a 'typical' null hypersurface, like a null-cone in Minkowski 
spacetime. In particular, this guarantees that, for some function k, defined on 
the horizon, £ a V a £* = n£ h . There will have much more to say about k further 
along. 

Now, consider the Null Raychaudhuri equation: 

£ a V a 6 {e) = -hf e) - a ab a ab + u ab u ab - 8ivT ab £ a £ b . (7.0.1) 

Where a ab = P a m P b n V \ m £ n) - \q a b0{i), and u ab = P a m P b n V [m £ n] . Since 
we know that the horizon is defined by the fact that dm = 0, and that it has 
£ a as a null tangent vector, we know that the left hand side of this equation 
vanishes on the horizon. Furthermore, we know that u ab vanishes, since the 
3 One might wonder about this statement regarding area. See (6.1) 
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pullback of £ a onto the horizon also vanishes. Thus, the above equation reduces 
to: 

= a ab a ab + 8itT ab £ a e b . (7.0.2) 

Since both of these terms are manifestly positive (the horizon metric is 
non-negative-definite, and the null energy condition is satisfied), this means 
that they both must independently vanish. Now, since the decomposition of 
P a m P h n VJ b = VJ b is: 



Va4 =^qab0(£) + (Tab + U ab , (7.0.3) 

=<t„6,. (7.0.4) 

and since i b V ' a i b = along with oo ab = a ab = 6(g) = 0, it must be the case that 
all contractions of q ma q n bV a £ b = V m £ n = 0. Therefore, we must have, on the 
horizon only, 

V a £ b = uJ b + u a k b . 

Since £ a has a fixed norm, however, we know that lcSb£ a = 0, since it is 
equal to minus itself by the product rule. Therefore, contracting the above 
expression for V a £ b onto £ b gives us the result that uj a = 0. Therefore, on the 
horizon, we have the condition: 

WJ b = u a £ b . . (7.0.5) 
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Where co a is called the Hajicek one-form, and should be taken as a 
three component object intrinsic to the horizon. Also, it should be clear that 
k = l a uj a - In computations, it is often simpler to take advantage of the fact 
that £ a k a = —1, and work with cu a = —k\S a t h since, among other things, the 



meaning of the pullback operator is clear, as one can compute the action of 
the derivative operator, perform the index contraction and then perform the 
pullback operation. 

Now, note that all of the above is completely invariant if we replace 
(£ a , k a ) with (L a , K a ) using the rescaling transformation: 



Since all results above were merely dependent on the inner product £ a k a 
being equal to minus one, and each of them being a null vector, and both of 
these are true of you rescale £ a by an arbitrary function f. Note, however, that 
when you apply this rescaling transformation, you induce a change in uj a : 



L a = f£, 



a 



Ka — j,k a ■ 



(7.0.6) 



K b V a L h = 



jhV a (f£ b ) , 

hV a £ b -\Vaf. 



(7.0.7) 




Which, after pulling back onto the horizon, gives us the result 



UJ a -> UJ a + VJn(f) . 



(7.0.9) 
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And, since k = £ a u a , we also get k — > fn + £ej ■ This rescaling freedom 
gives us the ability to, without loss of generality, impose more stringent con- 
ditions upon our horizon. In particular, if we choose an appropriate function, 
we can impose the condition that u a not evolve along flowlines of l a . More 
formally, if we make this choice, we can declare that the horizon be a "Weakly 
Isolated Horizon" (WIH) by saying that it satisfies: 

= [£ t , W a }& = £ e (u a £ b ) = £ b £ e u a . (7.0.10) 

Where we twice used the condition that £^£ a = 0. Now, the role of 
the rescaling freedom becomes clear. If u a does not satisfy (7.0.10), then we 
simply choose to rescale £ a and k a by some / according to (7.0.6), and solve the 
differential equation for / in such a way that (7.0.10) is satisfied. This will then 
fix the function / up to a positive constant (since l a needs to remain future- 
pointing and nonzero). Once this choice is made, it will impose the condition 
that k be constant on the Horizon (though it does not fix the nonzero value 
of k). This result is called the Zeroth Law of Black Hole Dynamics. 

7.1 The First Law of Black Hole Dyanmics 
7.1.1 A Simple Proof of the First Law 

We can use the isolated and dynamical horizon rules to derive the first law 
of black hole dynamics. The proof will be somewhat laborious, however, and 
rather than directly going into it, we will instead derive the first law quickly 
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here, and then generate the full proof. The below proof will have the advantage 
of being direct and easy. It will have the disadvantage of being very dependent 
on the details of the Kerr solution. This simpler proof is outlined in [15] 

Now, consider the Kerr spacetime, with metric as given in Appendix A. 
We define the horizon to be the surface at which A— 2 Mr = r 2 +a 2 — 2 M r = 0. 
Each closed 2-dimensional section will also have a constant t coordinate. Using 
the decomposition formalism built up in Chapter 6, it is apparent that the 2- 
dimensional metric of this surface is given by: 



/ A 2 + sin 2 {e) \ 



B 4 



(7.i. i; 



Where the coordinates are (0, 9) and A + and B + denote that the func- 
tions A and B as defined in Appendix A take the values of r required by the 
equation r 2 — 2 Mr — a 2 = 0. It is now easy to compute the area of the 
2-dimensional surface: 



a 2 



d 2 x ^/q= d9 d(f)A + sin (9) = 4vr A + , 
Jo Jo 

=4tt (r 2 + a 2 ) = An (V/ 2 + 2 MVM 2 - a 2 + M 2 - a 2 + 
=4tt (2 M 2 + 2 MVM 2 - a 2 ) . (7.1.2) 

Meanwhile, it is easy to calculate the angular momentum in the Kerr 
spacetime. Remembering the boundary formulation given in (5.2.2), we see 
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that it is necessary to add a term § d 2 x P l fii to the ADM Hamiltonian. Where 
P a is given by the solution valued integral: 




d 2 xf j (K ij -f j K) . 



(7.1.3) 



In order to find the angular momentum, we look for the <p component 
of this vector. Due to the vagaries of vector valued integrals, we must multiply 
by a unit vector in the <fi direction under the integral. After multiplying by 
this vector, and then taking the limit to infinity, we get: 

M y ' 

Now, we insert this relationship into the expression (7.1.2) to eliminate 
the dependence on a and replace it with dependence on J. The result is: 



Now, we take the variation of this equation, treating M, J, as indepen- 
dent variables with independent variations: 



1 



A = M 2 + VM 4 - J 2 . 



(7.1.5) 
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2MVM 4 - J 2 + 2M 3 



1 tA ( VM 2 - a? 
5 A ' 



vr V 2M 2 + 2MVM 2 - a 2 
a 



+ 5J 
1 



8tt 



2M 2 + 2MVM 2 -a 2 J ' 
«(JA + fi(JJ. (7.1.6) 



Where the multiplier k of ^-SA is the surface gravity of the horizon, 



8tt 

1 7-r7. j. _ j-i l : l nn.iJr}T~i v 

■b; 

while Q is the rotation parameter to the horizon, most directly defined as the 



defined from the two null normals l a and k a to the horizon by K = —£ a k b, V a i 



value that ^ takes on the horizon 1 . The generalization to the Kerr-Reissner- 
Nordstrom case is straightforward, and in the end will only involve the addition 
of a Q5Q term, where Q is the charge parameter of the black hole, and $ is 
the electrostatic potential at the horizon. 

Several things should be noted about this derivation. First, it was heav- 
ily dependent on the Kerr solution. We had to explicitly calculate the area of 
the Kerr horizon, and then take its variation. Therefore, this can really only be 
interpreted as a statement about the phase space of Kerr (Reissner-Nordstrom) 
spacetimes, and perturbations thereof, and not a general statement about 
black holes. Furthermore, it required some odd mixtures of terms-the mass, 
charge and angular momentum, in standard formulations, are properly defined 
by integrals at infinity. We performed one of these above in order to obtain 
the expression J = aM. On the other hand, Q, k, $, and A are proper- 



1 And, thus, in a na'ive sense, the "rate of rotation" of the space around the horizon. 
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ties inherently intrinsic to the horizon itself. They would, in principle, be 
computable even in spacetimes without boundary so long as those spacetimes 
contained closed trapped surfaces. It seems weird to have a single equation 
that multiplies these objects together haphazardly without much regard to 
which surface they belong. Older, more general formulations of the first law 
also have found themselves beset by this problem. And it is the task of the 
Isolated and Dynamical Horizon formulation to remove these problems by 
generating a formulation of the First Law that is both not dependent on the 
details of a particular solution and also in terms of quantities wholly defined 
on the horizon. In order to do this, however, we must continue further into 
the Isolated/Dynamical Horizon rabbit hole in order to build up some tools. 
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Chapter 8 

Decomposition of quantities on 
the horizon into ADM quantities 

Another key use of the IH/DH formalism is the ability that it gives you to 
decompose 3+1 quantities into quantities well defined on the horizon. Now, 
consider the case of either an isolated or dynamical horizon. In general, a 
spacelike section of a 3+1 split will intersect an isolated or dynamical horizon 
in 2-dimensional spacelike slices. Each of these slices will therefore have a 
spacelike normal s a lying in the 3+1 slice, and will also be perpendicular to 
the timelike normal n a used to perform the 3+1 slicing. Therefore, we can 
make the usual choice i a = a (n a + s a ) to decompose the outgoing null vector 
to the horizon, where the proportionality to a is chosen so that £ a mirrors 
the time vector t a as closely as possible. We will first develop some generic 
formalism that works for all zero expansion surfaces. Then, we will work out 
some quantities specialized to the Isolated and Dynamical cases. 
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8.1 Doubly decomposing ADM quantities on the 
horizon 

It should be clear that that the project described above, where one defines 
an ADM splitting, and then, within the context of that splitting, divides the 
3-space into the tangents and normals to an intersection of the horizon with 
the particular 3+1 slice, is exactly the sort of the double foliation discussed 
in Appendix D. The formalism developed above can now be used to derive 
powerful relationships between the intrinsic horizon quantities and the intrinsic 
3+1 quantities. In order to accomplish this, we are going to need to work out 
a few identities. First, recall from equation (3.1.16) that 

V n 6 = -n a V b ln (a) - K ah 

It is our goal to decompose the quantity V a s& in the same way. In order to 
do this, we write down the quantity, and then expand the 4-metric in terms 
of projection operators and normal operators: 

V a s 6 =g a c 9b d V c s d , 

= (la - n a n c ) (j b d - n b n d ) V c s d , 

=lalb d ^cSd - la C n b n d V c s d - n a n c % d V c s d + n a n b n c n d V c s d , 

=V a Sf> + la C n b s d V c n d - n a n b n c s d V c n d - n a n c ^ b d V c s d , 

=V a s b - n b s d K a d - n a n b s d V d ln (a) - n a j b d n c V c s d . (8.1.1) 
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So, we have decomposed the first three terms relatively easily in terms 
of 3+1 quantities. Furthermore, since the argument used in deriving equation 
(3.1.16) made no assumptions about the dimension of the space involved nor 
the signature of the normal vector, we can further decompose the first term 
into —H a b — s a Vbln , where H a b is the extrinsic curvature of the horizon 
section in the 3+1 slice, and B is the generalized Lapse function of the 2+1 
foliation in the neighborhood of the horizon 1 . The last term above, however, 
requires a more intricate analysis. First, we need to calculate the inner product 
of the unbarred radial vector with the normal vector 2 in terms of 3+1 variables: 



n's. =9"V»6 = aBg" b (V„r) (V b R) , 

= {is)' n ~ (8 - L2) 

Now, we compute the value of — 7fc d n c V ' c Sd- Before moving on, note 
that the fixed norm of s a means that the contraction of this term onto s b 
vanishes. Therefore, we can replace the 3-metric appearing in this term with 
the 2-metric. Now, we compute this term, taking advantage of the lack of 

1 Less formally, B = —r=^, where R =constant determines the location of the horizon, 

and is chosen as a coordinate on the 3+1 slice 

2 See the first section of this chapter for more clarity here. In essence, the unbarred radial 
vector is s a = BV a R, where B is chosen in such a way as to make s a a unit normal relative 
to g ab while s a — BV a R is a unit normal relative to 7 afc . 
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torsion to interchange the order of derivatives of a scalar: 



-lb d n c V c s d = - q bd n c V c (B^V e R) , 

= - q bd n c BY e V c V e R - q bd n c B (V e R) V cl de 

-q bd n c ^ (V e R) V C B, 
= - q b e n c BV c V e R - q bd n c B (V e R) (n d V c n e + n e V c n d ) 

-q bd n c s d V c ln(B) , 
= - q b e n c BW e W c R - B q bd n c n e (V e R) V ' c n d , 

= -B q b e W e (n c V c R) + B q b e ( V e n c ) V C R - B q bd (n e V e R) n c W c n d , 
= -B q b e V e (^n c s^j - q b c s d K cd - ^ q bd n c s c V d ln (a) , 

= - Bq b c V c ( -^^Sa) - q b c s d K cd - -f3 a s a V b ln (a) , 
\ aB J a 



= -Bq b 



V c « aa 



1 



q b c s d K cd -—B a s a V b a, 



a' 



=-f3 a s a V b ln (B) - - V 6 (f3 a s a ) - q b c s d K cd , 
a a 



1 



a 



V b ln (B) - -V b (3 ± - q b c s d K, 



a 



cd ■ 



(8.1.3) 



where we made the definition (3± = (3 a s a . Putting all of the above expressions 
together, we find that 
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V a s b = - H ab - s a V b ln (B) - n b s c K ac - n a n b s c V c ln (a) 

+ —n a V h ln (B) - -n a V b (3 ± - n a q b c s d K cd . (8.1.4) 
a a 

So, this allows us to completely decompose V a s b into ADM variables, 
while equation (3.1.16) enables us to decompose V a n b into ADM variables. 
We can now use this formalism to express the intrinsic horizon quantities in 
terms of ADM variables. Before fully doing this, we follow the notation of [2] 3 
and decompose the extrinsic curvature into components normal and tangential 
to the horizon: 



Jab =q a c qb d Kcd J a = q a b s c K bc K ± = s a s b K ab , 

K ab =la C lbK cd = {q a c + s a s c ) (q b d + s b s d ) K cd , 

=J ab + J a s b + s a J b + s a s b K ± . (8.1.5) 

Now, we are ready to work through the relevant intrinsic and extrinsic 

quantities. We start with the expansions, as they are the easiest. 

3 Though it should be noted that Cook does not define K± and instead just writes it as 
K — J. Also, note that Cook works under a gauge where B = (3± — s a f3 a . This is not 
consistent with a na'ive decomposition of the induced geometry on the horizon from the 
4-geometry, since if t is the time coordinate and r is the radial coordinate of the horizon, 
then B — y=?F and f3± — Ba 2 g tr , which are clearly two inequivalent expressions in the 
general case. 
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=q ab V a t b = q ah V a [a (n a + s a )\ , 
=a (q ab V a n b + q ab V a s b ) , 

=a(-J-H). (8.1.6) 

A nearly identical derivation will show that 0(k) = j^. {H — J). Imme- 
diately, we can therefore see that, in 3+1 language, the Isolated/Dynamical 
horizon condition is a requirement that H = —J and J > on the horizon. 
We can perform a similar decomposition to the shear: 



cr a b =q a c qb d ^Jd - ^q a bO , 

=q a c qb d V c [a (n d + s d )\ + q ab (J + H) , 
=a{-J ab -H ab )+ l -aq ab {J + H) . (8.1.7) 

This then shows us that the isolated horizon condition that a ab van- 
ish is equivalent to the condition that H ab = —J ab - Tracing this equation 
will then automatically give you the H = —J requirement for 9 = 0. Note 
that the twist u ab = q a c qb d ^ \J-d\ is automatically zero. This should not be 
surprising since the 2-dimensional horizon is smoothly embedded in both the 
3-surface and the 4-surface, and therefore must have a vanishing twist by 
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Frobenius's theorem 4 . Alternately, it is easy to substitute the definitions 
of n a and s a in terms of gradients of functions and then take both of their 
twists, and then show that the twists are proportional to V[ a V&]/ for some 
function in both cases, and are therefore zero. The Raychaudhuri equa- 
tion (6.2.2) is then satisfied if, in addition to J ab = —H ab , we also have 
T ab £ a £ b = a 2 (T ab n a n b + 2T ab n a s b + T ab s a s b ) = a 2 (p - 2j a s a + S ab s a s b ) = 0. 
These conditions are then sufficient to ensure that the horizon is a WIH [16] 
as defined in equation (7.0.10). 

8.1.1 Dynamical Horizon case 

Now, let us consider the case of a dynamical horizon (DH). These horizons are 

still defined as a topologically compact level surface of 9 = 0. If the horizon 

is expanding, however, it is necessarily the case that outgoing null vectors 

can only momentarily have zero expansion relative to infinity, because the 

expansion of the horizon will make the old horizon radius lie in the interior of 

the horizon at later times. In other words, the light ray hovering stationary 

on the horizon at time t will be falling inward toward r = at time t + St. 

Therefore, the outgoing null vector is no longer a tangent to the Dynamical 

Horizon, which goes from being a null surface to a spacelike surface 5 . By this 
4 See Appendix B.3 in [ 3], amongst other sources 

5 Contracting horizons will be timelike surfaces. They have several properties, however, 
that are not horizon-like. In particular, timelike 8 — surfaces will not be trapping surfaces- 
some null and timelike objects from inside these sorts of horizons will be able to escape to 
infinity. Furthermore, their existence necessarily requires the presence of matter violating 
the dominant energy condition. This may in fact be the case for black holes emitting 
Hawking radiation, but is otherwise not considered physically realistic. Therefore, we will 
not consider these surfaces here. 
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argument, then, the formerly null generators of the dynamical horizon must 
now be spacelike vectors. It is logical to adjust these generators in such a way 
as to make them as close to the i a as possible. We don't expect there to be 
dependence on the angular variables on the horizon, so the most logical choice 
for the generating vector v a is given by the ansatz: 

v a = t + ^s a . (8.1.8) 

It is easy to show that v a v a = 2 atp + tp 2 , which is positive so long as —2 a < ip. 
So, therefore, for positive ip, v a is a spacelike vector that corresponds to l a for 
■0 = 0. Now, we enforce the condition that v a generates a surface of 9 = 0: 

o = v c v c e = [i c + ^s c ) v c e , 

£ e e + ^s c V c (q ab V a £ b ) , 
s c V c [-a (J + H)\ , 

- s c ( J + H) V c a - as c V c J - as c V c H , 
s c eV c ln (a) - as c V c J - as c V c H . (8.1.9) 

Further simplification of the second and third terms on the right hand 
side above is involved enough that each term will be dealt with individually. 
Before we begin, recall the Gauss-Codazzi condition from equation (3.1.21) 
tells us, in the case of a 2-surface embedded in a 3-space, we have, after using 



£ v 6 



(£ v 9 — £t 
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the Hamiltonian constraint (3.1.23): 



R =R — 2s a s b R ab — H ab H ab + H 2 , 
2s a s b R ab —R — R — H ab H ab + H 2 , 

= (I6vrp + K ah K ab -K 2 )-R- H ab H ab + H 2 , 
= 16np + J ab J ab + 2J a J a + K\- (K± + Jf-R- H ab H ab + H 2 , 
s a s b R ab = X - (l67rp + J ab J ab - H ab H ab + 2J a J a - J 2 + H 2 - 2K ± J - . 

(8.1.10) 

Now, knowing this identity, we work toward decomposing s a V a H: 



s a V a H = - s c V c (q ab V a s b ) , 

= - S c q ab VSaS b ~ s c {V c q ab ) V a s b , 

= - s c q ab R cab d s d - s c q ab V a V c s b + s a s c (V a s b ) V a s b 

+ s c s b (V c s a ) V a s b , 
=s c s d R cd - q ab V a {s c V c s b ) + q ab (V a s c ) V c s b 
V b ln (B)] V b ln (§) , 

lQnp + J ab J ab - H ab H ab + 2J a J a -J 2 + H 2 - 2K ± J - R 



+ 

1 
2 



q ab V a 
1 



V b ln (B) 



H ab H ab + ^[Y"B)\ l ,B 



~-8np + ^ ( J ab J ab + H ab H ab + H 2 -J 2 -R)+J a J a - K ± J 
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(8.1.11) 



This takes care of the third term in equation (8.1.9). In order to take 
care of the second term, we first start with the momentum constraint, equation 
(3.1.26), and contract it onto s a 6 : 



8nfs a =s a V b (K ab - r b K) , 

=s a V b (J ab + s a J a + s b J a + s a s b K ± ) - s a V a (q bc K bc + s b s c K bc ) , 

= - J ab V b s a + V b J b + s a s b V b J a + K ± V b s b + s b V b K ± 
- s a V a J - s a V a K ± , 

=J ab H ab + V a J a + 2s a s b V a J b - K ± H - s a V a J , 
s a W a J =J ab H ab + V a J a - 2s a J b W a s b - K±H - 8nfs a , 

=J ab H ab + V a J a + 2J b V b ln (B) - K±H - M] a s a . (8.1.12) 

Therefore, putting equations (8.1.12) and (8.1.11) into equation (8.1.9), 

we get 



(£ v 9 - £ e 0) =s c 9V c ln (a) - a (j ab H ah + V a J a + 2 J b V b ln (B) - K±H 
-8nj a s a ) + Unp + l - (j ab J ab + H ab H ab + H 2 - J 2 - R 



6 Here, we make use of the decomposition in equation (8.1.5), and drop any terms that 
vanish either by the fact that two tensors are normal to each other or that vanish by virtue 
of s a having a fixed norm. 
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+J a J a -K ± J + ^V 2 B 



B 



--s c 9V r ln (a) — a 



(j ab + H ab ) (J ab + H ab ) -K ± (J + H) 



+ ^ (H + J)(H-J) + V a J a + J a J a - ±R + 8tt (p - fs a ) 

+2J b V b ln (B) + 4v 2 S 
B 



-s c 6VJn (a) - a 



- jfl + V„ J" + J"J„ + 8ir (p - s"j„) + 2J»V„(n (B) 

+Iv*b 



5.1.13) 



And finally, requiring that £ v vanish, and that the surface is a dynam- 
ical horizon where the outgoing expansion 9 vanishes, we get 



ip =- 



o 



2^ (v ab <Tab) ~\R + VaJ a + J a J a + 8tT (p - j a S a ) + 2 J a V Jn (B) + iV 2 £ 



-8TrT ab £ a e b - a ab a a 



a 



2^ (^ a6 ^) - |fl + V a J« + J a J a + 8n(p- j a s a ) + 2 J a VJn (B) + iV 2 5 

(8.1.14) 

Several things should be noted about the above equation. The first 
thing is that if) has a strong dependence on angular variations on the horizon. 
The traditional definition of the horizon angular momentum [11] [1] [17] is 
given by L = ^- § d 2 x ^/qK ab f a <f) a = ^ § d 2 x ^JqJ a <f) a . Secondly, note that 
the only derivatives of the 2+1 'lapse' function B appearing in (8.1.14) are 
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2-dimensional derivatives-apparently if) is not sensitive to radial dependencies 
in B 7 . 

Finally, one might be concerned by the presence of R in the denominator 
of this equation, since its value is heavily restricted by the topology of the 
horizon by the Gauss-Bonet theorem. But, in fact, the presence of R is actually 
quite fortunate, as in the case of the Schwarzschild metric, it is the only term in 
equation (8.1.14) that is not explicitly zero, and therefore, without the presence 
of R, we would have no control over ip approaching zero as our spacetime 
approached the spherically symmetric, non-dynamical case. 8 

8.1.2 Area Balance Law 

We now show how this formalism can generate an area balance law, and in 
turn, how this area balance law dictates conditions upon the lapse function 
and shift vector, as originally proved by Ashtekar and Krishnan [16]. The key 
insight in this derivation is to note that the Dynamical Horizon is a spacelike 3- 
surface locally defined by the level surface 9 = of the locally defined function 
9. Therefore, it is a well-defined Cauchy surface that can be chosen for the 
3+1 splitting of spacetime. Now, having made this definition, the next step 
is to calculate the mass flow across the horizon, and then to make what have 

now become commonplace 2+1 decompositions of the relevant geometrical 

7 This is something that we might expect, considering that B is defined as ,_ 1 =. 

Any radial variation in B is gauge, as it can be adjusted merely by a r paramaterization of 
the function R. 

8 For example, in the analysis of the Vaidya metric below, note that the M — > limit 
would be very poorly behaved without the 2-curvature term. 
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quantities. As we go through the following derivation, we will freely use the 
dynamical horizon condition 9 = to eliminate any terms of the form J + H 9 : 



d 3 x8irT ab £ a n b = J ^d 3 xa (MT ab n a n b + 8nT ab n a s b ) 



(R - K ab K ab + K 2 ) 



7 d 3 x a 



R + 2s a s R ab + H a H ab — H — J a J ab — 2J a J a 

-Kl + K\ + J 2 — 2K±J) - s a V b (J ab + J a s b + J b s a + s a s b K ± ) 
+s a V a (J + K A 



7 d 3 x a 



s a sY d Rcadb - K ± J - J a J a + ^{R + H ab H ab - J ab J ab 
-J ab H ab - s a s b V b J a - V b J b - s b V b K ± - K ± H + s a V a J + s a V a K ± ] 



^ d 3 x a 



s a V c V a s c - s a V a V c s c - J a J a - V a J a 



+ -[K + n ab H —J J ab — ZJ H ab 



2s a s b V b J a + s a V a J 



7 d 3 x a 



V c (s a W a s c ) + s c V a H + {V a s a f - {V a s c ) V c s a - J a J a 



-W a J a + =-(R + H ab H ab - J ab J at> - 2J a0 H ab ) + 2s° J a W b s 



ab 



rab : 



9 Note that the derivation to follow is done intrinsically to the Dynamical Horizon, using 
the constraint equations on the horizon. Therefore, it is best to think of the s a that appears 
below to be the unit normal parallel to the v a discussed in this section, and for the unit 
timelike normal to be the timelike normal perpendicular to this-namely, the unit normal of 
a timelike observer falling into the black hole with no velocity transverse to the black hole 
surface. 
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+s a V a J] , 

J ^d 3 xa [v a {s b V b s a ) - J a J a - V a J a - 2J a VJn (B) 
+ \{R- H ab H ab - J ab J ab - 2J ab H ab ^ + s a V a (J + H) 



(8.1.15) 



And now, we are nearly done. To further simplify this expression, we follow 
Ashtekar et al. [16], and make a few simplifying assumptions. First, we re- 
member that s a is tangent to the DH. Therefore, s a V a (J + H) = 0, and we 
can ignore the last term in (8.1.15). Second, we make a coordinate choice on 
the DH that makes B depend only upon the radial coordinate. This elimi- 
nates the total divergence term as well as V a B. Finally, we make a choice 
of the lapse function so that a = ~, which makes a^fy = y/q. This enables 
us to rewrite the 3-integral J d 3 x as J dr <f> ^fqd 2 x. The integral of V a J a be- 
comes an integral of a total divergence over a boundary, while the integral of 
R becomes the curvature invariant \ of the 2-dimensional section of the DH. 
Consequently, we can now rewrite (8.1.15) as 



J d 3 x8irT ab tn b = J dr j> ^dd 2 x -J a J a + ^ ab cr ab 

X J dr= J d 3 x (l6irT ab £ a n b + 2J a J a + a ab a ab ) . . (8.1.16) 

The right hand side of 8.1.16 is known as the Hawking energy, which has 
several interesting properties: 

First, if T ab satisfies an energy condition such that T ab £ a n b > 0, as one 
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would expect for a fluid with a positive density greater in magnitude than its 
pressure, then the right hand side of the equation is explicitly positive-definite, 
which then immediately tells us that r must be increasing and that % > 0, 
meaning that the apparent horizons that make up the DH must individually 
have spherical topologies, and that the DH satisfies the law of nondecreasing 
area as one would expect from the second law of black hole dynamics. 

Second, it is already apparent that (8.1.16) also contains the content 
of the first law of black hole dynamics. The left hand side term describes an 
increase in areal radius, which can be taken to be an infinitesimal area increase 
that can be identified with 5 A. The stress energy tensor term exactly describes 
the amount of mass entering the black hole, and thus can certainly be identified 
with 8M. Meanwhile, the last two terms involve nondiagonal terms of K a b, 
which are known to be associated with black hole angular momentum. There- 
fore, it shouldn't be very surprising that with some minor tweaks, (8.1.16) can 
be modified to produce a proof of the First Law of Black Hole dynamics, but 
with the advantage of being defined explicitly on the horizon, and that it is 
proved without the dependence on the details of the Kerr solution used in the 
proof given in (7.1.6). 

8.2 Dynamics of the Vaiyda solution 

Now, let us apply this formalism to one of the simplest dynamical examples: 
the Vaidya solution [18]. This solution describes a spherically symmetric null 
dust in the background of a Schwarzschild black hole. In order to make the re- 
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lationship between the 3+1 split and these boundary conditions more explicit, 
we are going to explicitly provide a 3+1 split, and then compute the necessary 
value for ip such that (8.1.9) is satisfied. 

Start with the Vaidya metric as it is commonly given, where M(v) is 
an arbitrary function of v: 



ds 2 = - M - 2A ^ ) J dv 2 + 2dvdr + r 2 d9 2 + r 2 sm 2 (6) dcf> 2 . (8.2.1) 

While this form gives us the simplest and most direct description of the 
metric, it is not ideally suited for a 3+1 split, as a computation of the inverse 
metric tensor will show us that dv a dv\,g ah = 0. Therefore, we transform this 
to Kerr-like coordinates by making the coordinate transformation v — t + r. 
Upon completing this transformation, the metric tensor takes the form (in 
(t,r,8,<f)) coordinates): 



/ _ ( i _ 2 



9ab 



2M(t+r) 



\ 



2M(t+r) 
r 






1 + 



2M(t+r) 
r 

















(8.2.2) 



r 2 sin 2 (9) j 



Where M is now an arbitrary function of t+r, rather than v. This metric 
has the exact same form as the Schwarzschild metric does in Kerr coordinates. 
A labourious but straightforward computation gives that the Ricci scalar of 
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this metric is equal to zero, and that: 



G a b — Rah — 



2M 



ab 



/ 1 1 ^ 
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V 



(8.2.3) 



/ 



Furthermore, we can compute two null vectors that have only t and r 
components, and after fixing their inner product with each other to be equal 
to —1, we find that they are: 



k a dx a =— (—dt — dr) , and 
2 c 



ell a dx a =c 



2M\ , / 2M 
-dt[l I + dr 1 



(8.2.4) 



where c is an as of now arbitrary constant. It is then easy to show that 
Qab = 9ab + Lh + kj b = r 2 dQ 2 , which then gives: 



q a b 



2M(t+r) 



2M(t+r) 
r 






^ | 2M(t+r) j q 



2M(U r- 






^ 



r 2 
r 2 sin 2 (9) J 
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/ ( 1 _ 2M(*+m 
2M(t+r) 



+ 



2 M(t+r) 



1 + 



V 



2 M(f +r; 
r 






r 2 d9 2 + r 2 sin 2 (0) # 2 . 













(8.2.5) 



After doing this, we then can easily find the trace against q ab the gradi- 
ent of any one-form that has only temporal and radial dependence, since the 
contraction on q ab will eliminate anything but angular terms, and the angular 
terms depend only on r: 



q ab V a v b =q ab d a v b -q ab T ab c v c , 

=0 - ^q ab v c (g ac ,b + 9bc,a - gab,c) , 

= ^q ab V C 9ab,c, 

= -v r . (8.2.6) 
r 

Now, we can use (8.2.6) to essentially read off the expansions of the two null 
vecors. The result is that 9(k) = ~y~ c an d = — (l — ^y L )- Since, at the 
point where r = 2M(t, r) we have 9^) negative and 9 zero, it is clear that this 
surface represents a trapped surface. Note that it is also clearly dynamical, 
due to the explicit inclusion of the time coordinate above. 

Now, we wish to 3+1 split this metric tensor. We will make the obvious 
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choice of t =constant slices and proceed from there. With this choice, we 
clearly have a = — 
tensor is given by: 



clearly have a = , 1 = = , (5 a dx a = ^y-dr, and n a = (—a, 0, 0, 0). The 3-metric 

V r 



( 1+ 2_M q 

r 



lab 



\ 



Or 2 
r 2 sin 2 (6) J 

while the extrinsic curvature can be calculated according to the rule 



V 



(8.2.7) 



Kab = ~ 7a C V c n fe , 

= - {lalb d d c n d - i d a lb e ?de c n c ) , 



Kn =0 - ry a . 



(8.2.8) 



which gives 




AM M 2M 6 M 2 



12.9) 



With all other extrinsic curvature components equal to zero. We now 
2+1 decompose onto the horizon. Our radial vector is given by 
r a = y^l + ^■(l, 0, 0), therefore giving us a spacelike "lapse" of ^/l + 
The 2-metric is simply the ordinary metric of a 2-sphere. Since K ab only 
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has diagonal terms, we get J a = 0, Kj_ = 1 \ M K rr and Jab 



2M 



-dQ 2 . 



Finally, we have the extrinsic curvature of the horizon in the 3-surface: 



H ab = ~ q a c qb d V c r d , 



Hab =^AB C r c 



: ^r C {^Ac,B + lBc,A ~ lAB,c) , 
11 ( d 



1 + 



2 M 



dr 



-.dQ 2 . 



1+ 2M 

1 ** 



(8.2.10) 



Therefore, on the horizon, we have Hab = —Jab, which tells us, au- 
tomatically, that a ab = and 8 = 0, the latter fact we, of course, already 
derived in the full 4- space. The vanishing of the expansion and of the shear 
tremendously simplifies (8.1.13), which, on the horizon, now becomes: 



(£ v 9 — £0) = —a 



-\ + 8n (p-ff a ) 



(8.2.11' 



A glimpse at the Raychaudhuri equation (6.2.2), coupled with all of 
these computations, will show one that the only term remaining in £$ is 



8nT ab £ a £ b . Setting £ v 9 = 0, and solving for ■0, 



if, 



8ttT, 



„6< 



l)a. fib 



a 



(i - 8ixT ab n a n b - MT ab n a f b ) ' 



(8.2.12) 
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We now remember that G a b = [dt + dr) 2 , which gives us G a bn a n b = 1 7Tu > 
G ab n a f b = -G ab n a n b and G afe ^ 6 = Putting all of this together into 

(8.2.12), we get the result 

(8.2.13) 

a 

Now, all that is left is to set the value of c by requiring that l a = 
a(n a + f ), so that £ a becomes the best approximation possible of the time 
evolution vector t a = an a + f3 a on the horizon. To do this, we simply calculate 
£ a — an a , and then require that the answer be orthogonal to n a 



n a (t-an a ) =£ a n b g ab + a , 



1 + 



2M 



1 + 



2M 



2M 



2M + 2M 



+ 



1 + 



2A1 



2c+ 1 



(8.2.14) 



So, after all of this, we find that the value ip takes on the horizon is 
4M^2. We can now use this quantity to calculate the time rate of change of 
the horizon's area: 

First, we calculate the tangent vector to the dynamical horizon: 
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v a =t + ^f a = c 
=d t + 4Md r . 



2M\ n ( 2M . 

l + —)dt+[l )O r 



+ 1p 



I _|_ 2M 

r 



<9 r , 
(8.2.15) 



where the hat indicates that the equality is only valid on the horizon. 

Now, it is easy enough to directly calculate the rate of area increase of 
the black hole: 



5 A —£ v j) d 2 x yfq , 

= jd 2 x l -^q ab v a d a q ab) 

= j d 2 xAM^r 2 sin9, 
=16ttM(2M) , 

=32ttMM. (8.2.16) 

Which is a logical conclusion, considering that the straightforward, di- 
rect computation of the black hole area gives us the value A = 4 irr 2 = 16 7rM 2 , 
whose first derivative is clearly 32 nMM . Clearly, this method was massive 
overkill for a situations such as the Vaiyda metric with a known analytical so- 
lution. However, in a dynamical spacetime, where the location of the horizon 
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is not known, this technique gives an exact way to trace not only the horizon, 
but also the normals to the horizon as it expands. 10 



10 Note that the area of the dynamical horizon is a slicing dependent quantity. See [19] 
for a direct analysis of the slicing dependence of the area of the Vaiyda DH. 
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Appendix A 

Null decomposition of the 
Minkowski and Kerr spacetimes 

In order to further clarify the above procedure, here, we will work out a cou- 
ple of examples of null decompositions. First, let us examine the Minkowski 
spacetime, given in spherical coordinates (t,r,9,4>): 

g ab dx a dx b = -dt 2 + dr 2 + r 2 d6 2 + r 2 sin 2 (6)d<p 2 . (A.0.1) 

We wish to examine the null geometry of the surface whose vector space 
is tangent to the outgoing vector £ a = (1, 1,0,0). To do this, we define v = 
—t + r, which gives g ab dv a = £ b . Then, we do a coordinate transformation, 
replacing t with v, which gives the metric: 
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9ab 



-110 

10 
r 2 

y r 2 sin 2 (0) J 



^010 
110 

oo^ 

I a ) 



(A.0.2) 



Where the first row/column represents a, and the next three coordinates 

are (r,Q,4>) as they normally would be. Obviously, da a = (1,0,0,0), and 

action upon this form with g ab gives £ a = (0, 1, 0, 0). We can see that, at any 

point in the spacetime, the tangent space is spanned by £ a , (0,0, -,0) and 
1 r 

(0,0,0, t-— ), along with one additional vector, normal to the second two 

r sin(0) 

vectors. We shall call this vector k a , and a simple computation shows that 
there is only one choice of k a that is null, normal to the two angular directions, 
and satisfies k a da a = —1. This choice is given by k a = ( — 1, — |, 0, 0), which, 
upon lowering with g a b, is equivalent to k a = (|, — 1, 0, 0) 

Now, we define the outgoing light cone at some time by setting v = 
constant, which means that dv = 0. Eliminating the appropriate columns 
from the above equation then gives the induced metric: 



/ 



Qab 



\ 



(A.0.3) 




Or 2 
^00 r 2 sin 2 (6) J 

Which is manifestly degenerate. The general matrix q ab satisfying q a mQ. mn (lnb 
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q ab is given by: 



q ab = 



\ 



a b 
c 



c 



1 



(A.0.4) 



2 (0) J 



But, we now require that q ab also satisfies q ab k a = 0. By inspection, we 

<— 

can see that k a = (— 1,0, 0). Therefore, we can see that q ab k a = — \{a,b, c), 

<- «- 

which then gives us: 



/ 



q ab = 







4 












\ 



(A.0.5) 



2 W / 



Which is, of course, the answer that one would have na'ively guessed 
to be the correct one. As before, one should carefully note that the the first 
index of the lowered induced metric acts upon what was the r index in the 
4-dimensional spacetime. As shown above, this means that the corresponding 
vector is £ a . Meanwhile, the raised indices also have an index that is labeled 
as r. This index, however, acts upon k a , and not upon £ a , which has zero 
pullback onto the null surface. Note (as stated above) that this is, in fact 
necessary, since the covector space to a vector space is defined precisely by 
the condition that the covector basis e , when acting on the vector basis e b 
must give the result 5 a b , which is impossible if £ a is a vector and if £ a is a 
covector. The simplicity of the above Minkowski example was chosen to show 
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the inevitability of this fact. 

Now, let us turn our attention the Kerr solution to the Einstein Equa- 
tion. This solution is an algebraically special solution, and therefore, can be 
given in the form of a fiducial Minkowski metric plus one of the principal null 
vectors times itself. In spheroidal coordinates (t,r,<f>,9), we have: 

9ab = Vab + Ck a k b g ab = rf - Ck a k b . (A.0.6) 



7- / -, B . . 2/1 „ x , / A Aasin 2 9 
ri ab = diag(-l,j,Asin 2 9,B) k a = (--,-l, , 0) . (A.0.7) 



C= 2M ^ B A = r 2 + a 2 B = r 2 + a 2 cos 2 9. (A.0.8) 



Where a and M are taken to be two parameters, with \a\ < M. It is 
easy to check that k a is null according to the Minkowski background metric, 
and therefore, it is also easy to verify that k a is null relative to the full metric 
as well. Therefore, either i] ab or g ab can be used to raise or lower indices on k a , 
and we take rf h to be the metric inverse of r\ ab . We are now going to investigate 
the properties of this spacetime's horizon, and the properties of its embedding 
into the full 4-dimensional spacetime. Before we do so, it is useful to introduce 
the other principal null vector of this spacetime: 
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A — 2Mr B(A + 2Mr) (A — 2Mr)asin 2 6 
a ~ (_( ^4 } ' 2A 2 ' 2A ' } ' 

Which, when raised using g ab gives: 



(A.0.9) 



A + 2Mr A — 2Mr {A + 2Mr)a 



(A.0.10) 



2A ' 2A ' 2A 2 
It can be verified that this vector has zero norm, and that £ a k a = £ a k a = 
— 1. Now, note that this vector seems to have some quite strange behaviour 
when A — 2Mr = 0. Noting from the definition above that A = r 2 + a 2 , it is 
easy enough to prove that r 2 — 2 M r + a 2 = is satisfied when: 



r ± = M± \[W 



(A.0.11) 



Note that, for these values of r, £ a oc dr, while £ a has a vanishing 
coordinate in the d r direction. We therefore expect that for these values of r, 
this is, in fact, a null submanifold with associated null vector i a . Using the 
above condition, and taking r = r + as the associated function (and thereby 
setting dr = 0), we take the 4- metric, do some algebra, and find the following 
induced metric: 



q a b 



a 2 sin 2 9 


— A+asin 2 9 





B+ 


B+ 


—A+asin 2 8 


A 2 + sin 2 9 





B+ 


B+ 











\ 



(A.0.12) 
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Where a subscript of + indicates that a function of r is taking on the 
value r = r + in its r argument. This is a more complicated looking metric 
than the one that we derived for the Minkowski spacetime, but, once again, 
it is easy to verify that it has vanishing determinant, is nonnegative definite, 
and that its zero eigenvector is (1, ^,0) = (£ t + , £^ + ,£ 9 +), and, since £ r + = 0, 
we can therefore just project this vector onto our null space and call it £ a . 
We therefore have a null tangent space to our horizon spanned by £ a and two 
spacelike angular directions. Finally, we pull back k a on the horizon with result 

7 _ (_A+ A + sin 2 e n \ 
— V B , B+ ' V r 

Solving the equation q a mQ mn Qnb = Qab gives the following matrix, for 
undetermined functions X, Y, Z: 



/ B++2A + a(sin 2 e)X-A 2 (sin 2 e)Y 
a 2 sin? 6 



V 



x 
z 



x z 



1 

B+ 



\ 



I 



(A.0.13) 



Contracting this onto k a and requiring that the answer be the zero 
vector then gives us the solution: 



X = 



Y = 



Z = 0. 



B+ B + sin 2 9 
Substituting this answer into the above expression for q ab gives us: 



(A.0.14) 
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q ab = 



\ 



a 2 sin 2 6 


a 


B+ 


B+ 


a 


1 


B+ 


B+sin 2 e 












1 

B~ + J 



(A.0.15) 



Which we can then take to be the inverse 3-metric on the Kerr Hori- 
zon. As a final note, it is often easier to work with the so-called "untwisted" 
coordinates on the horizon. These coordinates are defined by the coordinate 
transformation: 



at 



(A.0.16) 



Under this transformation, we can follow the above procedure, and 
obtain the following 3-metric (in (t,<p,0) coordinates): 



/ 



q a b = 











B A 



q ah = 



J 



( a 



\ 



a 2 sin 2 9 


a 





B+ 


A+ 


a 
A+ 


B+ 

A\sin 2 d 











1 

B+ 



(A.0.17) 

As one might guess, this transformation makes i a = (1, 0, 0) while leav- 
ing k a equal to the more complicated expression (—1, A +™ in e , 0). Having done 
all of this work, it is now simple enough to proceed using quantities defined 
intrinsically on the Kerr Horizon. Note that the above procedures are non- 
trivial in the Kerr case-the inverse of the (9, (j>) part of q a b is not equivalent 



to q 



ab 
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Now, having completed this Appendix, it should be clear how to go 
about inducing enveloping geometries onto null subspaces. In the end, the 
process is the same as for non-null subspaces, only with an extra step, because 
the tangent space and the cotangent space have to be treated separately. 
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Appendix B 

Hamiltonian of Klein- Gordon 
Field in an external gravitational 

field 

Here, I will decompose the action for a Klein-Gordon field coupled to an exter- 
nal gravitational field using a 3+1 formalism, and then derive the Hamiltonian 
and the associated equations of motion. In particular, we will show that while 
you can get boundary terms in the Hamiltonian, you get none of these bound- 
ary terms in the equations of motion. While it is not explicitly included here, 
the entire argument below is essentially unaltered if one were to add a self- 
interaction term V(<f>) to the action, so long as it dependent only on the value 
of the field and none of its derivatives. 

First, consider the Klein-Gordon action: 
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S = f d 4 x^\g\L = J d A x^\g\ (-Iv a 0V a - ^m 2 2 J . (B.0.1) 

Since we are concerning ourselves with only the dynamics of the scalar 
field in the background metric, and therefore holding the external metric fixed, 
and since we furthermore are not computing the stress-energy tensor, we will 
treat the square root of the determinant as part of the volume element, and 
treat the Lagrangian density as a true geometric density. This will simplify 
the calculation of boundary terms and does not change the content of any 
results. Therefore, from here on, the term involving the square root of the 
metric will be factored into the d A x term, which will be understood to have 
the appropriate weight to be a volume element. 

Knowing the Lagrangian density, it is easy to compute the momentum 
conjugate to the field: 

11 = = ^ (~l Va< ^ " 2 LmV ) = " VV • (B -°- 2) 
If this calculation were being done in Minkowski spacetime, then it 
would be direct to equate V'0 to — 0. Since we are instead in a 3+1 split, it 
is necessary to insert a factor of the inverse metric and decompose, yielding: 

-VV = - g ta V a ^ = - (-- 2 5t a + - 2 A V a , 
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= 10-1/3^0. 
cr cr 



(B.0.3) 



Where a is the lapse function and f3 l is the shift vector for the particular 
slicing in question, k... indicate 3- dimensional indices, and Vj denotes the 
connection compatible with the 3-dimensional metric 7^. Putting this all 
together, we get: 

n = ^ (0 - FVtf) . (B.0.4) 
So, having done this, now we can see that: 



= ^n0-^V a 0V l 0-^m 2 < /> 2 , 

= ^0 - ^/3>V,0 - \ (V - V.0V.0 - ^m 2 2 



ilty - ^ (/3*V,0) + ^ (/3 4 V,0) 2 - ^V>V,0 - ^m 2 2 



Now, we can combine (B.0.5) with the definition of the Hamiltonian 
density H = 110 — L to obtain the Hamiltonian density, and then we can use 
(B.0.2) to eliminate all dependence on 0: 
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X =110 - L 



2 ^ 2a 2 V ^ ; Y 2a 2 



1 (w) 2 + ivw + i m y 



« 2 - 1 - - 1 

:-n 2 + n (/3 j v^) + -v>v^ + -m 2 </> 2 . 



i 



2a 2 



(B.0.6) 



This then immediately lets us integrate the Hamiltonian density to get 
the Hamiltonian (with the factor of ypy once again an implied subcomponent 
of the d 3 x: 



H = / Hd 6 x 



d 3 x n 2 + n (/? l v,0) + l - V>V,,0 + ^m 2 ^ . 

(B.0.7) 

Which gives the equivalent phase space action 



S 



d x 



2 j. 2 

m (p 



(B.0.8) 



And we therefore take (B.0.7) to define our Hamiltonian for this Klein- 
Gordon system. It is now trivial to get the equation of motion for </>: 



127 



d x 



110- 

fl + « 2 n + p'v. 



y n 2 + n (^v<0) + \ v>v,0 + ^m 2 2 



(B.0.9) 



Which, since the variation must be equal to zero under an arbitrary 
variation, gives us the equation of motion <fi = a 2 U + f3 l Vi<f). The variation 
with respect to required to obtain the equation of motion for IT is somewhat 
more intricate however: 



5 S = ! r/'.r - 11 4 (urn - V V 



/^-n+(, 
+ ^ J d 4 x (Vi (n/3» - ^riVi/f - 0/^11) , 
o + / An = / <?x (mV - v<v,0 - - /?<v ( n) + / A (,n,-) , 
^/A[„V-v-v^-v,(n^)] + /A(^), 

(B.0.10) 

= J d 4 x(m 2 <f)-V i V i <i)) . (B.0.11) 

Where, in going from (B.0.10) to (B.0.11) we noted that the divergence 
in the bulk integral precisely canceled the boundary integral via a simple ap- 
plication of Gauss's theorem. Therefore, we obtain the equation of motion 
for II, which is simply n = V*Vj — m 2 0. We have now found the complete 
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equations of motion for this Hamiltonian system after taking a variation of 
our dynamical variables, and importantly, we have found that there are no 
dangling boundary terms remaining after our variation has been taken, mean- 
ing that it is not necessary to add any counterterms to our original action-the 
variation of the phase space action will be zero for any (II, 0) that satisfy the 
equations of motion for the Klein-Gordon field. 

As a final step, it would be remiss to not check our result against the 
traditional empty-space Klein-Gordon equations of motion in Cartesian coor- 
dinates. This is simply equivalent to setting 7^ = 8^, and making the choices 
— and a = 1. If we do this we have: 

= 11 IT = <9*<9;0 - m 2 • (B.0.12) 

If we take a time derivative of the equation for 0, and then substitute 
this answer into the equation for II, we get 

=d % di<p — m 2 (f) , 

= - + d l d t (f) - m 2 (f) , 

=7] ab d a d b (j) - m 2 (j) . (B.0.13) 
And (B.0.13), indeed, is the Klein-Gordon equation of motion. 
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Appendix C 

Boundary terms in the 
S chwar z schild- Rob er son- Walker 

spacetime 



As another example, let us consider the case of the Schwarzschild-Robertson- 
Walker spacetime. This is furnished by a simple extension of both the Schwarzschild 
and Robertson- Walker spacetimes, and is given by the line element 1 



2M\ 

g ab dx a dx b = - 1 )dt 2 + (a(t)Y 



2M 



+r 2 d6 2 + r 2 sin (6) d(j) 2 ] . (C.O.I) 



1 Note: the function a(t) is meant to be evocative of the Robertson- Walker function a{t). 
It is not meant to have anything to do with the Kerr parameter a. 



130 



If we take the limit a(t) = 1, a(t) — > d(t) — > we clearly recover the 
Schwarzschild spacetime, while if we take the limit M — > we recover the flat 
Robertson- Walker model. Furthermore, direct computation of the curvature 
generated by this line element shows that its Einstein tensor is: 

~ ^ (frfj) " ^sin> (0) . (C.0.2) 

Which manifestly satisfies G/ = Gg e = G/ = - (^^)- We can 
therefore interpret the line element given in equation (C.0.1) as the gravita- 
tional field of a mass M surrounded by a strain-free dynamical fluid with a 
radial current 2 . Note that this solution, however, is likely unphysical, as com- 
puting the trace of the Einstein tensor above will show that both the Ricci 
scalar and R a i ) R ab contain a singularity at what would otherwise be the black 
hole horizon at r = 2M . Assuming a polytropic equation of state for the fluid 
can produce a concrete form for the function a(t), and given this form, numeri- 
cal computation of the geodesic equations generated from equation (C.0.1) will 
show that this coordinate is accessible in a finite amount of proper time, so 
the black hole horizon in the McVittie solution is a naked singularity. Further- 

2 by changing the time coordinate to r = 3ln(a) + In (l — we can diagonize the 

Einstein tensor. This will make the fluid pressure in the r direction different from the 
fluid pressure in the angular directions as well as making the metric tensor nondiagonal. 
Therefore, the fluid is not isotropic in its comoving frame 



131 



more, computation of the scalar C a b c< iC abcd out of the metric (C.0.1) produces 
a result that is also singular, making it unlikely that this singularity arises 
solely due to a caustic in the matter 3 . 

Furthermore, varying the line element with respect to t will generate 
the equation: 



dEschwarzchild • ( -2 ^ , o2 2 i 12 2 • 2 / q\ \ 

-aa [ r + " r + 4 ) 7 sin (9) J , 



dt 



r 



E R bertson-W alker ■ (CO. 3) 

a 



Here Eschwarzschild is the formerly conserved quantity in the pure Schwarzschild 
spacetime arising from the former killing vector dt, and E Robertson-Walker is the 
conserved energy in the pure Robertson- Walker spacetime that is derived by 
combining the geodesic equation for the time coordinate with the unit timelike 
condition. 

In particular, this means that for an expanding (contracting) universe, 
geodesies will tend to have their energy redshifted (blueshifted) as time pro- 
gresses. This dynamical effect will make stable circular orbits impossible, and 
on a timescale governed by the Hubble time and by the mass of the central 
object, initially circular orbits will tend to fall into the singularity at r = 2M 
(or spiral outward toward infinity and become unbound). Numerical solutions 



3 There is one exception to this: if a(t) satisfies ^| = — which then makes a(t) an 

exponential function. In this case, C a b c dC abcd and R become finite, although R a ^R ah retains 
its singularity. So, while the singularity remains, its strength is lessened and its character 
changed. 
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of the geodesic equations verify this property of geodesies and generalize it to 
elliptical orbits as well. A few plots of solutions with polytropic equations of 
state are plotted below. The geodesic equation was solved numerically using 
the Mathematica software package, which also produced the plots. 

20 - 




-20 - 

Figure C.I: A graph of the solution for an equatorial Geodesic with Tq = 15M, 
Ho = 5,41 ^ f 10 — , and L = L C i rcu i ar . The red circle represents r = 2M, while the 
gold circle represents the innermost stable circular orbit at r = 6M. The test 
particle completes 21 full orbits before falling into the black hole. 

Since a particular orbit is only bound for a finite amount of time in this 
spacetime, we therefore expect to find some interesting time-dependence in the 
energetics of this spacetime. In particular, this example was chosen as a way 
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Figure C.2: A graph of the solution for an equatorial Geodesic with r = 15M, 
Hq = 5,41 ^ f 10 — , and L = l.lL circu i ar . An addition of a relatively small amount 
of additional angular momentum than this produces a simple escape orbit, 
approximately at L — 1.35 M. The red circle represents r = 2M, while the 
gold circle represents the innermost stable circular orbit at r = 6M. The test 
particle completes 21 full orbits before falling into the black hole. The test 
particle completes 28.9 full orbits before plunge happens 

to show that it is, in fact possible to generate reasonable-seeming spacetimes 
that have time dependent ADM masses, as shall be shown below. 

So now, let us compute this spacetime's ADM mass. For the slicing co- 
ordinate, we choose the coordinate t from the line element above. Meanwhile, 
for the fiducially flat comparison background metric 7^, we make the obvious 
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Figure C.3: A graph of the solution for an equatorial Geodesic with r = 15M, 
Hq = 5,41 ^ J L0 — , and L = L circu i ar . First, note that the Hubble constant for 
this plot is ten times the Hubble constant for the other two plots. The red 
circle represents r = 2M, while the gold circle represents the innermost stable 
circular orbit at r = 6M. The test particle completes 3 full orbits before falling 
into the black hole. Note the degree to which there is already a poor division 
between the inspiral phase and the plunge phase, although the geodesic does 
manage to complete three full orbits. Orbits with values of Hq much larger 
than this tend to not even complete whole orbits before plunge, even when 
given large amounts of angular momentum. 



choice of the spherical Euclidean metric given by choosing a(t) = 1, M = 1 . 

4 One might object to this choice, arguing that the spacetime is not, in fact, asymptotically 
flat. However, each timeslice under this choice IS asymptotically flat, and with a rescaling 
of the time coordinate, this spacetime is conformal to Schwarzschild spacetime, which is 
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Then, denoting covariant differentiation relative to 7^ as Vj, we compute the 
ADM mass: 



IQ-kMadm = 




And then, when we take the limit r — > 00, we obtain the result Madm = 
Ma 2 , indicating that the ADM mass of this spacetime is time-dependent. 

asymptotically flat. Also, since none of the derivatives in the ADM formula is a time 
derivative, and there are two factors of the lowered metric (inside the parentheses, and in 
the two-metric determinant), and one factor of the raised metric, this choice yields the same 
formula as you would obtain by choosing a fiducial metric by defining M = 
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Furthermore, this behaviour explains the odd behaviour of the geodesies of 
this spacetime-rather than interpreting the energy of an orbiting particle as 
being redshifted by the cosmological expansion, we can also interpret the mass 
concentrated in the central object increasing by an amount proportional to 
the scale factor over time. Therefore, initially circular orbital paths eventually 
find themselves attracted by a mass larger than the one they were attracted 
by initially, and see the radius of their orbit shrink. Eventually, they find 
themselves beyond the stable limit of their orbit, and have no choice but to 
plunge into the surface of the singularity. 

The conclusion to reach is that while, for asymptotically flat space- 
times, the ADM mass is conserved, there do exist examples where one can 
find a time-dependence in the ADM quantities. Relativity's insistence on lo- 
cal mass-energy conservation does not necessarily imply global mass-energy 
conservation. This result arises from the fact that this spacetime does not 
admit a global timelike Killing vector on the "sphere at infinity" 5 . 



5 In classical mechanics, it is not wholly abnormal for boundary conditions to change the 
overall value of the Hamiltonian, even if the bulk states stay the same. Consider a system 
of a string vibrating in its fundamental mode on a string of length L, mass per length fj, and 
amplitude A. Now, a simple evaluation of the Hamiltonian of this system [2( ] will show you 
that the total energy in the string is given by — • Therefore, if one were to lengthen the 
string without changing its tension, and if L were much smaller than the characteristic wave 
velocity in the string, one would expect that the string would stay in its new fundamental 
mode, and have energy • Thus the slow expansion of the string would allow the string 

to do work, causing it to lose energy. The proposal in this appendix is that we can similarly 
interpret time-dependent boundary terms in relativity-they are work terms indicating net 
global energy generated by the gravitational field, just as the factor expanding the classical 
string creates/removes net energy from that system. 
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Appendix D 

On various foliations of a spacelike 
surface in a 3+1 split 

Note that, when dealing with boundary conditions in 3+1 foliations, we often 
have to deal with TWO foliations of boundary terms-one involving the choice 
a constant r slice, giving a spacelike surface and one involving the choice of 
a constant R slice, which gives the boundary of that spacelike surface. In 
this appendix we will work out some consistency problems involving these two 
choices. The general scope of this argument is following that given in [ ], but 
the notation chosen in this derivation is chosen so as to match the notation 
in the rest of this work, as certain results derived in this appendix will be 
used heavily elsewhere in this work, particularly in Chapters 3, 5, and 6. 
Also, certain details are worked out in a slightly different way. For the below 
derivation, unless otherwise noted, we will be using the inclusion operator to 
map intrinsic 2- and 3-geometries into the enveloping 4-geometry, since we will 
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be comparing geometric quantities to each other. Also, we will freely raise and 
lower vectors using the full 4-metric g a b and its inverse. 

Now, consider a 4-manifold with boundary M with topology R x m 
with boundary m| Um|^.U(lRx dm). We will label the R portion of the 
spacetime with the function r. Similarly, there is some function R such that 
R is constant on dm and such that R is not constant on some neighborhood 
of dm. For now, we will assume that V a -R is a spacelike object, and examine 
the null case at the end of this appendix. Now, as was done in (2.1), we can 
easily define the two normals to this surface as they live in the enveloping 
4-space. And from these normals, we can construct projection operators onto 
their respective normal spaces: 

n a = «V fl r 

1 

a = — 

vV 6 (V a r) V fe r 

7 ab = g ab + n a n b 

Note, however, that it is not the case that n a s a = 0, since we have 
done nothing to guarantee that these things are normal to each other. In 
particular, if we choose r and R as coordinates, so that our coordinate system 
is (r, R, x 3 , x 4 ), it should be clear that it should not generally be the case that 
g rR = on any particular surface. Therefore, n a and s a are not a particularly 
good choice of basis for the normal vector space to dm. 



s a = BV a R, 



B 



vV fc (vjrj v b R 



(d.o.i; 



= ab = g ab _ s a s b 
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So, what to do? Clearly, it will be necessary to use the projection 
operators in a Gramm-Schmidt procedure. Since we are starting with the 
same 4-manifold and we are projecting down to the same final 2-manifold, 
it shouldn't matter which vector we start with and which vector we project. 
As we will see below, however, there is not complete equivalence 1 . So, we 
can define two different projections of the vectors, first onto the barred space 
defined by ^ ab and second, onto the double barred space defined by 7° . Finally, 
we define the projection operators defined by these two different methods, and 
note that they must be the same, as the target 2-space is simply <9m 



s a = BV a R = B~f a b V b R n a = a V a r = a l a b V b r , 

B=-2== a= 1 (D.0.2) 

vr b s a s b y/\r b n a n b \ 

g ab = g ab + n a n b _ -a-b = g ab _ ^6 + =a = b 

Now, the obvious conclusion that we can get from equations (D.0.2) is 
that n a n b — s a s b = —s a s b +n a n b . We are going to use this fact in order to derive 
an interesting conclusion regarding the two sets of vectors. In what follows, 



-"^For an extreme example for why this should be the case, consider the Schwarzschild 
spacetime in Kerr-Schild coordinates at r = 2M. If we slice by Kerr Schild time first, then 
we get a spacelike 3-manifold that stretches all of the way to the singularity at r = 0. The 
r = 2M surface appears like a perfectly well-behaved spacelike sphere. Null geometry would 
never have to be used or invoked. 

Now, however, if we were to choose to slice the spacelike slice FIRST, we would get a null 
3-surface, and only recover the spacelike 2-geometry after setting the null parameter of the 
horizon equal to a constant. Care must be taken in order to ensure consistency between 
these two approaches 
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note that the constructions of all of our projection operators are made in such 
a way that they annihilate relevant vectors. Therefore, n a s a = n a s a = 0. it is 
not the case, however, that the other contractions of these normals are going 
to vanish. Therefore, we will make the definition n a s a = ip. 
Expanding the definition of s a : 



s a = B la b V b R = B (5 a b + n a n b ) V h R , 
= B V a R + Bn a n b V b R , 

= g i S a + n a i\)) , 
B _ 

Sa= r^Sa-lpHa. (D.0.3) 



Similarly, 



n a = ^ a b V b T = a (5 a b - s a s b ) V fe r , 

= <5 (V a r - s a s b V b r) , 

= - (n a - s a ip) , 
a 

= -n a -ip- (Srs a - ipn a ) , 
a a \B j 

= | (1 + - (|f) - (D.0.4) 

So, we have now solved for the reverse order ADM vectors s a and n a 
in terms of their standard order counterparts, the "lapse" functions and the 
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parameter ip. Next, we enforce the fact that s a and n a is an orthonormal 
basis. Then, we will check to make sure that the consistency relationship 
derived from equating the 

q ab ig 

valid. Starting with equations (D.0.4) and 
(D.0.3), and assuming that all of the 'lapse functions' are positive 2 : 



s a s a =1 = (5jS a - ^n a j (^s a - im a ) , 
B 1 



B ^l + W' 
Continuing this process, 



(D.0.5) 



=\ 2 /=D\ 2 



, (i + ^ 2 ) + Fi ^ 

a ' \ aB 1 



= \ 2 



5 _ 1 _ S 
a + ip 2 ~ B 



[l + 2^ 2 + ^ 4 - (1 + ^ 2 )^ 2 ] , 

(D.0.6) 



While the condition n a s a = gives no new information. 
Now, substituting equation (D.0.6) back into equations (D.0.3) and 
(D.0.4) yields 



2 And if they are not, minus signs can be absorbed into the R and r in order to enforce 
this positivity 
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n a = ^l + ^ 2 n a - tps a , 
s a = y/l + ^ 2 s a - ^n a . 



(D.0.7) 



Which is clearly a Lorentz boost with boost parameter given by tfj = 
sinh(4>). Therefore, inverting the order in which one splits in a 1+1+2 split 
formalism will be equivalent to a boost when done the other way. 

And, as a consistency check, we compute: 



fi a fi b - s a s b = (yi + ^n a - (Vl +^ 2 n b - ^s 6 ) , 

= (1 + i) 2 ) n a n b - i>^l + i; 2 (n a s b + n b s a ) + ip 2 s a s b 

- (1 + V 2 ) s a s b + ip^l + ip 2 (s a n b + s b n a ) - i) 2 n a n b , 
=n a n b -s a s b . (D.0.8) 

So, as promised, our two projection operators are identical. 
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D.l Example: calculating the velocity parame- 
ter of Kerr Spacetime 

Now, for the sake of completeness, let us calculate the value of the boost 
parameter ip for the Kerr spacetime. Take the metric to be given by: 



9ab 



2 Mr 



2Mi 



2 Mr 



1 + 



B 

2 Mr 
B 



a sin 2 (9) 



(l+ 2 ^r)asm 2 (9) 




2 -^asin 2 {6) 
- {I + 2 -^) a sin 2 {6) 
A+ 2Mra l smHe) )stn 2 (6) 



B 



(d.i.i; 



,ab 



\ 





2Mr 
B 








2 Mr 
B 


A—2M r 
B 


a 
B 








a 
B 


1 

Bsin 2 (9) 














B 



(D.1.2) 



/ 



In coordinates labeled by (t,r,(j),9), and all other parameters labeled 
as in Appendix A. Take the function r = — t, and the function R = r. Then, 
we obtain via a simple reading from the above expression of the inverse to the 
metric tensor: 
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— - — ,0,0,0 



s ° =|0 ^zdLi7? ' '' 



ij) =n a s a = g ab n a s b 



2Mr 



B 



B 



B 



2Mr 



v /(fi + 2Mr) (A-2M>. 



B + 2Mr J V A-2Mr ' 

sinh{(f>). (D.1.3) 



Then, solving 1 = cosh 2 (<f)) — sinh 2 (<j)) gives: 



cosh(4>) 



A B + 2 M r a 2 sin 2 (6) 



(D.1.4) 



(A — 2 Mr) (B + 2Mr) ' 

Which immediately allows us, in analogy with Special Relativity, to find 
the velocity parameter v, 



v = tanh(<p) 



2Mr 



sinh(4>) 

cosh{4>) ~ ^/AB + 2Mra 2 sin 2 (6) ' 



(D.1.5) 




Figure D.l: The velocity parameter v plotted against r/M and 9. The plane 
is the surface / = 1 and a = 0.9. 



This term can be interpreted as the velocity of the 3-surface of constant 
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R relative to the 3-surface of constant r.. Note that it is zero on the surface 
at spacelike infinity, as r — > oo. Since g rt = ^ — > for asymptotically flat 
spacetimes, this result is only sensitive to asymptotic flatness, and not the 
details of the Kerr solution. Furthemore, for the case of asymptotically flat 
spacetimes, one can use the fact that v — at infinity in order to interpret v 
as the velocity of a point's reference frame with respect to conformal spacelike 
infinity. In order to bolster this interpretation, note that v takes the value 
of 1 for all points on either horizon r = M ± \JM 2 — a 2 , is greater than one 
in the region between the two horizons, and less than 1 everywhere else 3 . A 
Mathematica plot of v for a = .9 M in units of r/M is included in figure (D.l). 



3 As another aside, note that v takes the value for r = 0. Observers inside the inner 
horizon cannot access the region outside the horizon, making it somewhat nonsensical to 
talk about "local speeds of spacetime" relative to conformal infinity. One could, however, 
interpret v as the velocity of the reference frames relative to the point at the center of the 
Kerr ring. 
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